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I.  Summary 

An  investigation  was  undertaken  to  determine  if  by  cryomilling  NiAJ  and  TiB2  powders 
ultra  fine  subgrain  or  grain  size  material  could  be  obtain.  If  this  ultra  fine  grain  size 
composite  were  produced  then  it  was  assumed  there  should  be  a  significant  increase  in 
the  strength.  Cryomilling  did  decrease  the  grain  or  subgrain  size  but  the  strength  of  the 
20  V%  TiBj/NiAl  composite  was  the  same  as  20  V%  TiB2  made  by  the  XD™  process. 
Also,  a  0  V%  TiBj/NiAl  composite  produced  from  cryomilled  NiAl  was  about  the  same 
strength  as  the  20  V%  TiB^NiAl  composites.  The  above  mentioned  composites  were 
made  by  HIPPing  and  extrusion.  If  the  0  V%  TiBj/NiAl  cryomilled  composite  is  tested  in 
the  as  HIPPed  condition,  the  strength  is  much  greater. 

A  TEM  investigation  has  shown  that  in  case  of  extruded  cryomilled  composites  there 
is  a  dispersion  of  small  10-20  nm  Al203  and  AIN  particles.  In  the  case  of  the  as  HIPPed 
cryomilled  composite,  the  Al203  and  AIN  formed  a  network  at  the  grain  boundaries. 

It  may  be  speculated  that  the  greater  strength  of  the  as  HIPPed  composite  could  be 
due  to  the  network  of  Al203  and  AIN  acting  as  a  continuous  filament. 

The  various  theories  that  have  been  proposed  to  account  for  fracture  toughness  of 
discontinuous  metal  matrix  composites,  or  the  lack  there  of  ,  were  examined.  It  was 
observed  that  the  fracture  toughness  is  independent  of  particle  size  and  the  fracture 
proceeds  through  the  matrix.  A  significant  factor  which  has  not  been  considered  in  the 
past  is  that  the  matrix  is  in  a  severely  cold  worked  condition.  The  cold  worked  state  is  the 
major  factor  contributing  to  the  low  fracture  toughness  of  discontinuous  metal  matrix 
composites. 

The  Finite  Element  Method  (FEM)  was  employed  to  investigate  the  matrix  plastic  flow 
in  a  whisker  reinforced  SiC/AI  composite  under  external  tensile  load.  It  was  found  that  the 
plastic  zone  induced  by  the  plastic  relaxation  of  thermal  stresses  expands  under  the 
external  tensile  load.  The  overall  matrix  plastic  flow  was  characterized  by  the  expansion 
and  interconnection  of  the  plastic  zones  around  whiskers.  This  process  can  be  divided 
into  several  characteristic  stages,  and  related  to  the  global  stress-strain  relationship.  It 
was  also  found  that  composite  asymmetric  constitutive  response  to  external  tensile  and 
compressive  loads  as  well  as  the  composite  asymmetric  Bauschinger  behavior  was  due 
to  the  thermally  induced  plasticity  and  residual  stresses. 
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Interfaces  in  XD  Processed  TiB2/NiAI  Composites 


L.  WANG  and  R.J  ARSENAULT 

Interfaces  of  TiB:-NiAl  and  a-ALO,-NiAI  in  TiB:/NiAl  composites  have  been  investigated  by 
analytical  electron  microscopy.  Although  no  consistent  crystallographic  orientation  relationships 
have  been  found  between  NiAl  and  TiB:  or  ALO,.  semicoherent  interfaces  between  a-ALO, 
and  NiAl  have  been  observed  by  high-resolution  electron  microscopy  (HREM)  in  areas  where 
the  low  indexed  crystallographic  planes  of  a-ALO,  aligned  with  that  of  NiAl.  No  semicoherent 
interfaces  between  NiAl  and  TiB:  have  been  observed.  Silicon  segregation  was  consistently 
detected  by  X-ray  energy-dispersive  spectroscopy  (EDS)  at  the  TiB:/NiAl  interface  region. 
Segregation  has  not  been  detected  in  the  a-ALO,-NiAl  interface  region.  The  segregation  layer 
observed  at  the  TiB,-NiAl  interface  is  too  thin  to  absorb  any  of  the  thermal  residual  stress. 


I.  INTRODUCTION 

It  is  generally  accepted  that  excellent  bonding  be¬ 
tween  the  metal  matrix  and  the  ceramic  reinforcement 
can  be  achieved  by  the  XD*  (exothermic  dispersion) 

*XD  is  a  trademark  of  Martin  Marietta  Corporation.  Bethesda.  MI) 

process  for  several  of  the  metal  matrix  composites  rein¬ 
forced  with  TiB:  particles,  such  as  TiB.-NiAL1' :| 
TiB^-TiAL1'1  and  TiB.-Al  alloy.1’1  During  the  XD  pro¬ 
cess.  the  TiB:  particles  are  formed  in  situ  in  the  molten 
metal  matrix,  so  they  have  clean,  nonoxidi/ed  interfaces 
and  typically  are  single  crystals  of  high  purity.  Usually, 
the  crystal  structures  of  the  metal  matrix  and  the  ceramic 
particles  are  so  different  that  no  simple  topotaxial  rela¬ 
tionship  exists1'1  in  XD-processed  metal-matrix  compos¬ 
ites.  Recently,  it  was  found  that  during  high-temperature 
deformation,  dislocations  can  be  generated  from  a  matrix- 
particle  interface  into  both  the  matrix  and  the  particles 
in  XD  processed  TiB:/NiAl  composi'e.1' h|  Therefore, 
interface  structures  might  be  very  important  in  control¬ 
ling  the  mechanical  properties  at  these  composites. 
However,  another  intriguing  observation  concerning  the 
TiB:/NiAl  composites  is  the  lack  of  dislocation  gener¬ 
ation  due  to  cooling  from  the  high-temperature  ( 1673  K) 
annealing.1''’1  The  predicted  thermal  residual  stress  should 
be  very  large,  and  this  stress  should  be  relieved  by  dis¬ 
location  generation.171  It  is  possible  that  an  interface  layer 
exists  between  the  TiB:  and  the  NiAl  which  can  act  as 
a  stress  absorber.171 

The  purpose  of  the  current  investigation  was  to  deter¬ 
mine  if  an  interface  layer  exists  and  then  determine  if  it 
is  of  sufficient  thickness  and  also  to  determine  if  any 
consistent  crystallographic  relationships  could  be  found 
between  the  matrix  and  the  reinforcement. 

II.  EXPERIMENTAL  PROCEDURE 

Composites  of  0,  10.  and  20  vol  pet  TiB:/NiAI  were 
purchased  from  the  Martin  Marietta  Corporation 


L  WANCi.  (irailuiiic  Student,  and  R  J  ARSI'NADI.T.  Professor, 
are  with  the  Metallurgical  Materials  Laboratory.  Department  of 
Materials  and  Nuclear  Knginccring.  University  of  Maryland.  C  ollege 
Park.  MI)  20742-21  15 
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(Bethesda.  MD).  A  series  of  compression  tests  at  1033. 
1143.  and  1273  K  were  performed  in  the  strain-rate  range 
from  10  "to  10  :s  1 .  The  transmission  electron  micros¬ 
copy  (TEM)  foils  were  obtained  from  the  compressive- 
tested  sample  by  electric  discharge  machining  disks  w  ith 
a  thickness  of  0.5  mm  and  then  dimple  grounding  to  a 
thickness  of  0.07  to  0. 1  mm  at  the  center.  Foils  for  con¬ 
ventional  TEM  were  prepared  by  twin-jet  electrolytic 
polishing  in  a  5  pet  perchloric  acid-ethanol  solution, 
whereas  samples  for  high-resolution  electron  microscopy 
(HREM)  were  prepared  in  the  same  way  as  for  conven¬ 
tional  TEM  except  they  were  ion  milled  from  10  minutes 
to  1  hour  to  clean  the  surface  as  the  final  procedure. 
Most  of  the  conventional  TEM  work  was  carried  out  on 
the  l  MeV  high-voltage  electron  microscope  (HYEMi  at 
the  Argonne  National  Laboratory  (Argonne.  ID.  High- 
rcsolution  electron  microscopy  and  X-ray  energy- 
dispersive  spectroscopy  (EDS)  analyses  were  performed 
on  a  JEOL  2(XX)  FX  analytical  microscope.  The  X-ray 
wavelength  dispersive  spectroscopy  (WDS)  analyses  were 
carried  out  on  a  JEOL  X40A  scanning  electron  microscope. 

It  is  necessary  to  tilt  the  sample  to  such  an  orientation 
that  strong  low-index  reflections  exist  on  both  sides  of 
the  uueitace  and.  for  best  results,  the  electron  beam  is 
parallel  to  the  interface  when  performing  HREM  inves¬ 
tigations  of  these  composites  Although  there  is  a  very 
small  chance  of  obtaining  two-dimensional  (2-D)  lattice 
fringes  on  both  sides  of  the  interlace,  when  the  orien¬ 
tation  relationship  between  the  particles  and  the  matrix 
is  random,  it  is  usually  not  difficult  to  get  one-dimensional 
(l-D)  lattice  fringes  in  both  the  matrix  and  the  particles. 


III.  EXPERIMENTAL 
RESULTS  AND  DISCUSSION 

Figures  I  and  2  are  TEM  micrographs  of  0  and 
20  vol  pet  TiB  /NiAl  composites,  respectively,  which 
were  taken  on  the  1  MeV  HVEM  at  a  low  magnification, 
representing  the  general  view  of  the  microstructures  of 
the  composites.  Particles  which  appeared  in  the  0  vol  pet 
sample  arc  a-ALO,.  There  may  be  as  much  as  5  to 
10  vol  pet  of  these  particles  in  all  of  the  0.  10.  and  20  vol 
pet  TiB:/NiAI  composites.  Therefore,  in  Figure  2.  which 
is  a  micrograph  of  a  20  vol  pet  TiB-/!\’iAI  composite. 
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Hy.  i  TLM  micrograph  of  ()  vol  pet  TiB:/NiAI  composite.  Lou 
magnification. 


Hg.  2  TLM  micrograph  of  20  vol  pet  TiB:/NiAI  composite.  Lou 
magnification. 


there  are  particles  of  both  TiB;  and  a-Al20,.  Most  of  the 
larger  ones  are  TiB2  particles,  since  the  average  particle 
size  ofa-Al,0,  is  smaller  than  that  of  TiB:  particles.  The 
grain  size  and  particle  size,  including  that  of  a-ATO,. 
remain  constant  throughout  the  whole  deformation  pro¬ 
cess.  independent  of  the  testing  temperature  and  strain 
rate. 

It  has  been  found  that  the  matrix  composition  deviated 
from  the  stoichiometric  composition  by  quantitative  X-ray 
WDS  investigations,  i.e..  the  matrix  contains  —52  at. 
pet  nickel.  Apparently,  this  is  partly  due  to  the  formation 
of  the  a-AI:0,  particles  during  the  processing  which 
consumed  a  certain  amount  of  aluminum. 

An  investigation  was  undertaken  to  determine  if  any 
crystallographic  orientation  relationships  existed  be¬ 
tween  the  reinforcement  (including  «-AI;0,)  and  matrix 
materials.  Although,  in  some  cases,  one  or  two  lower 
indexed  crystal  planes  of  a-AI:0,  or  TiB:  phases  can  be 
aligned  with  that  of  the  matrix  phase,  no  consistent 
crystallographic  orientation  relationships  were  found  be¬ 
tween  the  reinforcements  and  the  matrix  phase.  In  other 
words,  there  was  a  random  crystallographic  orientation 
relationship  between  the  particles  and  the  matrix. 

However,  despite  the  differences  between  their  cry  stal 
structures  and  their  random  crystallographic  orientation 
relationships,  there  are  some  interesting  aspects  of  inter¬ 
faces  between  Al.O,  and  NiAl.  In  some  cases,  it  seems 
that  one  or  two  disordered  atomic  layers  on  the  a-AI:Q( 
side  exist,  as  shown  in  Figure  3.  It  is  possible  that  this 
observation  is  an  image  effect.  However,  it  has  been 
proposed  that  a  disordered  region  occurs  to  obtain  atomi¬ 
cally  well-matched  interfaces  where,  crystallo- 
graphically.  there  is  a  large  level  of  atomic  mismatch. |M 
Unfortunately,  there  are  at  least  two  difficulties  pre¬ 
venting  us  from  carrying  out  computer  image  simula¬ 
tions  to  interpret  the  observations.  The  first  difficulty  is 
that  the  crystal  orientations  between  NiAl  and  «-AI  ()< 
are  random.  As  a  consequence  of  this,  there  is  a  second 
difficulty:  i.e..  there  is  a  lack  of  theoretical  models  for 
interface  structures.  There  are  too  many  possible  way  s 
to  construct  the  interface.  Since  the  observed  disordered 
layer  is  apparently  in  «-AI:0.  which  has  a  higher  elastic 
modulus  than  NiAl.i:  “ 1,11  it  makes  the  interpretation  even 


Fiji  '  HKF.M  micrograph  of  <i-AI  O.-NiAl  interface.  Arrows  indicate  ihe  disturbed  interlace 
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more  difficult,  except  we  believe  that  elastic  anoma¬ 
lies1"  i;i  of  a-Al.O,  occur  at  the  interface  region. 

In  the  case  of  the  lower  indexed  planes  of  u-AEO, 
aligned  with  that  of  the  matrix,  semicoherent  interfaces 
were  observed.  For  instance.  Figure  4  showed  a  semi¬ 
coherent  interface  w  hen  the  {1(H)}  planes  of  NiAl  aligned 
with  {1 10}  (in  a  rhombohedral  coordinate  system)  planes 
of  a-AEO,.  By  using  a  hard  sphere  model  of  the  inter¬ 
faces  structure,  the  distance  between  misfit  dislocations 
is  [)  =  </,(/,/(</,  -  (/O.1"1  where  </,  and  </.  are  the  inter- 
planar  distances  for  reinforcements  and  the  matrix,  re¬ 
spectively.  The  interplanar  distances  of  ( 100)  plane  in 
NiAl  and  (110)  planes  in  «-AI:(),  are  0.2XX  and 
0.3479  nm.  respectively.  Calculated.  D  -  1.673  nm. 
which  is  just  about  six  times  as  large  as  the  interplanar 
distance  of  ( 1(H))  planes  in  NiAl.  In  Figure  4.  six  of  the 
( 1(H))  planes  in  NiAl  planes  aligned  with  five  of  the  ( 1 10) 
planes  in  a-AEOt  planes,  with  one  misfit  dislocation  on 
the  matrix  side.  This  agrees  with  the  calculated  value 
very  well.  Coherent  or  semicoherent  interfaces  between 
TiB;  and  NiAl  were  not  observed  via  HREM.  At  the 
TiB:-NiAI  interface,  there  exists  a  very'  thin  layer  (maybe 
a  different  crystal  structure  but  one  which  is  too  thin  to 
obtain  a  microdiffraction  pattern)  of  atoms  which  had  a 


different  lattice  spacing,  as  shown  in  Figure  5.  The 
thickness  of  this  layer  was  about  0.6  nm  when  the  ob¬ 
servation  was  made  along  the  basal  plane  of  TiB:. 

The  results  discussed  above  indicate  that  the  interface 
layer  at  the  TiB:-NiAI  interface  is  very  thin  (0.6  nm): 
therefore,  it  is  not  thick  enough  to  absorb  the  thermal 
residual  stress,  but  there  is  a  possibility  that  the  thermal 
residual  stress  is  relaxed  by  diffusion  of  atoms  which 
have  a  larger  atomic  radius  than  the  matrix  atoms  in  the 
interface  region.  This  relaxation  would  occur  because  the 
matrix  is  in  hydrostatic  tension1'1  and  the  oversized  atoms 
would  reduce  the  stress.  The  possibility  of  segregation 
was  investigated  by  X-ray  EDS  analysis.  In  Figures  6 
and  7.  which  are  TEM  images  of  a-AI.O«  and  TiB:  par¬ 
ticles.  respectively,  locations  A  through  C  are  the  areas 
at  which  the  EDS  micrographs  were  taken,  which  rep¬ 
resent  within  the  particles,  the  interface  region,  and  the 
matrix,  respectively.  A  beam  size  of  —20  nm  was  em¬ 
ployed  for  the  interface  region  and  a  slightly  larger  size 
for  the  particle  and  matrix  regions.  In  the  EDS  profiles, 
there  is  a  small  carbon  peak  which  is  due  to  contamination. 

Figure  5(a)  is  the  EDS  profile  taken  inside  an  u-AI:Ot 
particle.  Apparently,  there  are  no  impurities  inside  the 
o-Al-O,  particles  except  tor  a  trivial  amount  of  nickel 


|;ig.  4  HR  KM  micrograph  of  «- AM),- NiAl  interface  Arrows  indicate  the  misfit  dislocations  m  the  NiAl. 


I  ig  5  HRIM  micrograph  of  TiB  NiAl  interlace  Arrows  indicate  the  silicon  segregated  la\er 
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Fig.  6 — TKM  micrograph  of  «-AKO<  particles  of  0  vol  pet  TiB;/NiAl 
composite. 


which  is  not  much  more  than  background.  Profiles  taken 
in  the  AI:0,-NiAl  interface  region  and  the  NiAl  matrix 
close  to  the  interface  are  also  shown  in  Figures  8(b)  and 
(c).  A  small  amount  of  sulfur  was  detected  both  at  the 
interface  region  and  in  the  matrix.  In  comparing  the  counts 
obtained  from  the  interface  region  and  matrix,  it  is  sug¬ 
gested  that  there  is  no  preferential  segment  for  sulfur 
atoms.  Figure  9(a)  presents  the  profile  obtained  inside  a 
TiB2  particle.  TiB:  particles  are  as  pure  as  the  a-AEO, 
particles.  Figures  9(b)  and  (c)  are  profiles  taken  of  the 
TiB;-NiAI  interface  region  and  of  the  NiAl  matrix  close 
to  the  interface.  A  remarkable  silicon  peak  was  found  at 
the  interface  region  when  comparing  the  profile  obtained 
at  the  TiBj-NiAl  interface  with  those  obtained  from  the 
NiAl  matrix  or  TiB3  particles.  Since  a  silicon  peak  was 
only  detected  at  the  TiBrNiAl  interface  region  by  EDS, 
it  might  be  speculated  that  this  thin  layer,  as  observed 
in  the  HREM  investigation,  might  be  mainly  composed 


t  ig  7  —  TKM  micrograph  of  TiB;  particles  of  20  vol  pet  TiB./NiAl 
composite. 


(a) 


( b ) 


U) 

Fig  X — EDS  profiles  of  a  Al  t),  particles  and  close-by  regions: 
(<;)  inside  the  o-AI-O.  particle.  </>)  interface  region  of  ri-Al.O,  and 
NiAl.  and  <«-)  NiAl  matrix  close  to  the  o-Al.O,  panicle. 

of  silicon  atoms.  The  exact  structure  and  composition  of 
this  segregation  layer  are  unknown,  and  we  encountered 
the  same  difficulties  in  attempting  to  carry  out  a  com¬ 
puter  image  simulation  as  in  the  case  of  the  a-AI.O,-NiAI 
interface.  However,  the  consistency  of  the  observation 
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( b ) 


<c) 

f  ig.  V'-IiDS  profiles  of  TiB;  particles  ar.J  close  bv  regions: 
(</)  inside  ihe  TiB:  particles.  (/>)  interlace  n  ;ion  of  TiB  and  NiAl 
(silicon  A',,  peaks  arc  visible),  anil  (<  )  NiAl  matrix  close  to  the  TiB- 
particles 


of  the  interface  layer  and  Si  segregation  in  all  of  the 
samples  investigated  confirmed  the  existence  of  a  Si-rich 
interface  layer  at  the  TiB.-NiAl  interface.  Because  sili¬ 
con  is  one  of  the  most  common  impurities  in  commercial 
aluminum  metal,  the  silicon  atoms  segregated  at  TiB, -NiAl 


interfaces  would  probably  come  from  impurities  within 
the  aluminum  powders  which  are  used  in  the  processing 
of  the  composites.  Silicon  was  not  found  at  the  NiAl 
subgrain  or  grain  boundaries.  Obviously,  the  TiB_,-NiAl 
interfaces  are  energetically  preferential  sites  for  silicon 
atoms.  Experimentally,  coherent  or  semicoherent  inter¬ 
faces  between  NiAl  and  TiB,  have  not  been  observed 
via  HREM.  Generally,  interfaces  in  in  situ  composites 
have  low  energies  and  are  usually  semicoherent. |M| 
However,  the  tendency  for  forming  semicoherent  inter¬ 
faces  might  be  reduced  by  introducing  silicon  atoms  to 
the  TiB:-NiAl  interface  region. 

It  is  obvious  from  a  consideration  of  the  above  data 
that  an  interface  layer  does  not  exist  at  the  ALO.-NiAl 
interface.  The  Si  segregation  layer  at  the  TiB, -NiAl 
interface  is  too  thin  to  act  as  a  shock  absorber  of  the 
thermal  residual  stress. 

The  bond  strength  of  a  metal /ceramic  interlace  is  mainly 
determined  by  both  the  imperfections,  such  as  voids  and 
microcracks,  at  or  close  to  the  interface  and  the  work  of 
adhesion.115"'1  This  is  defined  as  W{li/  =  y,„  i-  y,  —  ym, 
where  y,„  and  y,  are  surface  energies  of  free  surfaces  of 
metal  and  ceramic,  respectively,  and  ym,  is  the  energy 
of  a  relaxed  metal/ceramic  interlace.  There  is  no  direct 
correlation  between  bond  strength  and  work  adhesion  in 
a  number  of  metal /ceramic  interlace  systems.  The  bond 
strength  seems  to  depend  mostly  on  the  interface  im¬ 
perfections  and  local  stresses.1151  For  the  TiB, /NiAl 
composite  samples  which  were  deformed  30  pet  at  a 
temperature  range  of  760  °C  to  I  (XX)  °C.  no  voids,  cracks 
or  any  other  debonding  evidence  has  been  observed  at 
the  interface  region,  and  the  HREM  investigations  have 
shown  that  both  a-Al:0<  and  TiB-  formed  a  perfect  atomic 
level  match  with  the  matrix  materials  (NiAl)  after  de¬ 
formation.  Therefore,  both  n-AIO,  and  TiB,  formed  an 
excellent  bond  with  the  matrix  (NiAl). 

IV.  CONCLUSIONS 

Based  on  the  above  observations,  the  following  con¬ 
clusions  can  be  reached: 

1 .  In  XD-processed  TiB,/NiAI  composites.  NiAl  formed 
an  excellent  bond  with  the  reinforcement  TiB,  par¬ 
ticles.  as  well  as  with  the  «-AI,0,  particles. 

2.  No  consistent  crystallographic  orientation  relation¬ 
ships  have  been  found  between  the  NiAl  matrix  and 
TiB;  or  «-ALO,  particles,  although  in  some  cases, 
one  or  two  low-indexed  crystallographic  planes  of  TiB, 
or  ar-ALO,  might  be  aligned  with  those  of  NiAl. 

3.  Silicon  segregation  has  been  found  in  the  TiB -NiAl 
interlace  region.  No  coherent  or  semicoherent  inter¬ 
faces  have  been  observed  between  the  TiB,  and  the 
NiAl.  and  the  Si  layer  is  too  thin  to  absorb  anv  of 
the  thermal  residual  stress. 

4.  There  is  a  strong  tendency  for  the  formation  of  semi¬ 
coherent  interfaces  between  o-AI-O,  and  NiAl.  No 
segregation  has  been  found  at  these  interfaces. 
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In  a  recent  paper.  Miller  and  Humphreys  (MH)[1]  make  several  statements 
which  should  be  examined  in  further  detail,  and  this  examination  will 
follow.  However,  there  i-  one  point  that  MH  make  which  should  be 
emphasized;  "However  as  the  ratio  of  fibre  length  to  diameter  decreases, 
the  approximations  made  in  the  modified  shear  lag  approach  become 
questionable. " 


First,  let  us  begin  with  the  general  concept  of  strengthening  of 
aluminum  and  aluminum  alloys  by  the  addition  of  whisker  or  particulate  Sic. 
MH  state  that  several  factors  should  be  included  when  considering 
the  strengthening  of  SiC/Al  composites.  Arsenault  et  al.  have  produced 
numerous  publications  of  the  list  of  the  factors  which  can  contribute  to 
strengthening  and  have  examined  the  contribution  to  strengthening  of  each 
one  of  them  in  detail  [2-31],  and  have  concluded  that  the  strengthening 
was  due  to  the  work  hardening  of  the  matrix.  For  example,  in  one  of  the 
reference  books  cited  by  MH,  there  is  such  a  list  given  by  Arsenault  [19]. 


Second,  let  us  now  consider  some  specific  points  raised  by  MH. 

One:  They  state  that  the  equation  derived  for  the  strengthening  due  to 
an  increase  in  dislocation  density  generated  by  the  difference  in  the 
coefficients  of  thermal  expansion  should  indicate  that  the  strengthening  is 
directly  related  to  volume  fraction  of  reinforcement  and  inversely  related 
to  particle  size  to  the  1/2  power.  Let  us  consider  the  equation  generated 
by  Arsenault  and  Shi  [8]  to  predict  the  increase  in  strength  due  to 
dislocation  generation  by  the  relaxation  of  the  thermal  residual  stress 
which  is  caused  by  the  difference  in  thermal  coefficients  of  expansion 
between  the  particulate  Sic  and  the  matrix.  This  particular  equation 
predicts  (as  stated  when  first  derived)  the  most  efficient  method  of  stress 
relaxation  due  to  dislocation  generation.  In  other  words,  this  equation 
predicts  the  increase  in  minimum  dislocation  density.  However,  more 
realistic  attempts  to  obtain  the  upper  limit  of  the  dislocation  density  due 
to  thermal  residual  stress  relaxation  have  resulted  in  an  unrealistically 
high  dislocation  density. 


Further,  in  regard  to  the  simple  equation  derived  by  Arsenault  and  Shi 
[8],  it  does  show  that  the  strength  change  due  to  dislocation  generation  is 
a  function  of  the  volume  fraction  of  the  particulate  divided  by  the  diame¬ 
ter  of  the  particulate  to  the  !j  power,  i.e.. 


Vf 

ho  «  - 


d* 


(1) 


Two:  MH  presume  that  dislocation  generation  will  only  occur  during  the 
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quenching  of  the  sample.  In  other  words,  if  the  sample  is  furnace-cooled 
(cooled  slower  than  l*C/min),  then  there  will  be  no  dislocation  generation. 
In  Fig.  1  (dislocation  density  vs.  volume  fraction)  and  Fig.  2  (dislocation 
density  vs.  particle  size)  the  data  was  obtained  from  samples  which  were 
FORNACB-COOLSD . 


0  5  10  15  20 

VOLUME  FRACTION  (%) 


Fig.  1  The  change  in 
dislocation  density  as  a 
function  of  volume 
fraction  of  SiC  (SiC 
whiskers  have  a  length 
to  diameter  ratio  on 
average  of  -  2  and  the 
dia  -  0.5  *im).  (Ref. 
27). 


PARTICLE  SIZE  (D) 


Fig.  2  The  change  in 
dislocation  density  as  a 
function  of  particulate 
size.  (Ref.  2/.) 
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If  we  consider  the  data  in  Figs.  1  and  2  and  make  the  assumption  that 
the  dislocation  density  is  uniform  (which  it  is  not)  and  use  the  following 
equation 

Ao  =  apbd.p'/'  (2) 

where  a  is  a  constant  ranging  from  0.5  to  1.25,  u  is  the  shear  modulus,  b 
is  the  Burgers  vector,  and  A p  is  the  change  in  dislocation  density  due  to 
the  ACTE  effect,  then  it  is  possible  the  predicted  increase  in  strength  is 
due  to  the  increase  in  dislocation  density.  The  predicted  increases  in 
strength  are  given  in  Tables  I  and  II. 

Three:  MH  state  that  the  change  in  subgrain  size  due  Sic  addition  to 

A1  is  inconsequential.  Figures  1  and  2  are  experimental  measurements  of 
the  decrease  in  subgrain  size  with  increasing  volume  fraction  and  an 
increase  in  subgrain  size  with  increasing  particle  size.  Employing  this 
experimental  data  with  the  experimental  data  of  McQueen  et  al.  [32]  who 
measured  the  change  in  strength  of  Al  as  a  function  of  subgrain  size,  it  is 
possible  to  determine  the  strength  increase  due  to  the  decrease  in  subgrain 
size.  As  can  be  seen  in  Tables  I  and  II,  there  can  be  a  considerable 
increase  in  the  strength,  for  a  particle  size  of  0.5  ym,  AoSg  =  131  MPa. 

Four:  MH  state  that  the  increase  stress  i.e.  work  hardening  between 

the  proportional  limit  and  the  yield  stress  defined  at  0.2%  offset  strain 
is  due  to  the  formation  of  Orowon  loops  about  the  Sic  whiskers  or 
particles.  It  has  been  shown  by  a  detailed  FEM  analysis  that  the  plastic 
zones  around  the  reinforcement  expand  upon  loading [33].  There  is  no 
indication  of  deformation  beginning  in  the  matrix  and  expanding  towards  the 
reinforcement. 

Fifth:  Let  us  now  consider  the  following  statement  of  MH,  "The 
strength  of  particulate  MMCs  is  consistent  with  the  estimates  based  on 
micromechanical  models  developed  by  conventional  2-phase  alloys."  The 
question  that  follows:  Which  conventional  model  has  taken  into  account  the 
strengthening  due  to  dislocation  generation  and  the  reduction  in  subgrain 
size  due  to  the  ACTE  effect?  As  far  as  we  know,  the  models  put  forward  by 
Arsenault  et  al.  were  the  first  to  take  into  account  the  strengthening  due 
to  an  increase  in  dislocation  density  CAUSED  BY  THE  ACTE  EFFECT  and  the 
reduction  in  subgrain  size  DOE  TO  THE  ACTE  EFFECT. 

In  conclusion,  several  points  need  to  be  made: 

•  The  strength  increase  due  to  the  addition  of  SiC  whiskers  or  particu¬ 
late  to  aluminum  or  aluminum  alloys  can  more  than  adequately  be 
accounted  for  by  the  increase  in  dislocation  density  and  the  reduction 
of  the  subgrain  size  due  to  ACTE  effect. 


•  Arsenault  et  al.  have  considered  several  other  possible  contributing 
factors  and  have  shown  that  they  are  inconsequential.  In  the  past, 
Arsenault  et  al.  have  proposed  that  the  strength  of  SiC/Al  composites 
is  due  to  a  nonconventional  phenomenon,  i.e.,  the  increase  in  disloca¬ 
tion  density  and  the  decrease  in  the  subgrain  size  due  to  the  ACTE 
effect. 
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Table  z 

Predicted  Increase  in  Yield  Strength  of  20  VV  sic  1100  Matrix 
and  Experimentally  Measured  Values 


Particle  Size  t\o  p 

D  (tm)  (MPa) 

AoSG 

(MPa) 

aores 

(MPa) 

A°ypred 

(MPa) 

Aov 

yexp 

(MPa) 

0.5 

126.6 

131 

0 

257 . 6 

153 

9 

99 

69 

17 

151 

88 

70 

79 

20 

17 

82 

69 

250 

59.7 

0-13 . 8 

17 

42.7-56.5 

26 

Aop 

Tha  predicted 

increase 

due 

to  the 

increased  disloca- 

tion  density 

AoSG 

The  predicted 

increase 

due 

to  the 

reduced  subgrain 

size 

^rt.s 

The  predicted 

average 

tensile  residual  stress 

^^pred 

The  sum  of  &o p 

+  A(7SG 

-  Aores 

A°yexp 

The  yield  stress  defined  at 

0.2%  : 

plastic  strain 

minus 

the  0.2%  yield 

stress 

of  the  0  V% 

1100  matrix 

[Ref.  27] 


Table  II 


Predicted  Increase 
Function  of 

in  Yield  Stress  of 
Volume  Fraction  of 

6061/Al  Alloy  as  a 
Sic  Whisker 

Volume 

AoP 

AoSG 

Acres 

Aoypred  Aoyexp 

V%  MPa 

MPa 

MPa 

MPa 

MPa 

1 

37.7 

0 

1.7 

36.0  6.9 

5 

56.6 

7 . 0 

8.6 

55.0  52 

20 

70 

48.3 

34.5 

83.8  100 

hap 

The  predicted 
tion  density 

increase 

due  to  the 

increased  disloca- 

floSG 

The  predicted 
size 

increase 

due  to  the 

reduced  subgrain 

Aores 

The  predicted 

average 

tensile  residual  stress 

AoYpred 

The  sum  of  6op 

+  Aosg 

“  Aares 

A°yexp 

The  yield  stress  defined  at  0.2% 

plastic  strain  minus 

the  0.2%  yield 

stress 

of  the  0  v% 

1100  matrix 

[Ref.  27] 
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ABSTRACT:  An  analytical  investigation  ot  the  dtllerenee  in  the  con¬ 
stitutive  behavior  of  silicon  carbide-aluminum  (SiC-AI)  composites 
under  uniaxial  compressive  and  tensile  loading  was  earned  out.  and 
the  associated  changes  of  the  residual  stresses  were  studied  It  is 
suggested  that  the  observed  asymmetric  response  ot  the  constitutive 
behavior  is  primarily  due  to  the  existence  of  residual  stresses  in  the 
composite.  These  residual  stresses  are  initially  introduced  because  ot 
the  differences  in  the  coefficients  ot  thermal  expansion  I  AC  1  I  i  dur 
mg  the  cooling  process  and  subsequently  changed  alter  the  external 
load  is  applied.  A  two-dimensional  (2-D)  finite  element  analvsis  ot 
a  hexagonal  array  ol  Si(  whiskers  in  an  Al  matrix  was  performed, 
the  result  of  which  testifies  to  the  appropriateness  ot  the  proposed 
explanation  It  is  concluded  that  the  2-1)  finite  element  analysis  is  an 
economical  way  of  adequately  reproducing  the  most  prominent  lea- 
lures  of  the  material  s  constitutive  behavior.  Bawd  upon  the  theory 
of  mechanics  ol  composite  material,  simplified  analytical  models  were 
developed  which  can  be  used  to  investigate  the  influences  ol  the 
whisker-matrix  interface  shear  load  transfer  and  the  volume  mismatch 
on  the  residual  stresses  It  is  concluded  that  the  mechanism  ol  load 
transfer  in  terms  of  normal  stress  at  the  whisker  lip.  which  i'  governed 
by  the  volume  (bulk)  mismatch  between  the  matrix  and  the  whisker, 
is  predominantly  responsible  for  altering  the  specific  pattern  ol  the 
residual  stresses  under  the  applied  load  In  contrast,  the  cited  ol 
shear  load  transfer  at  the  longitudinal  whisker-matrix  interlace  is  only 
restricted  to  a  very  small  region  in  the  vicinity  of  the  whisker  and 
does  not  influence  the  overall  symmetry  ol  the  constitutive  behavior 

KKY  WORDS:  thermal  residual  stress,  deformation  induced  residual 
stress,  asymmetric  constitutive  behavior,  whisker-reinlorced  metal- 
matrix  composites.  FEM.  volume  mismatch,  interface  load  translcr 
(shear  lag) 

Nomenclature 

tr„  Stress 

(cr„>  Average  stress 
A*  Area  of  element  K 
r?„  Effective  average  stress 
q  Applied  load  intensity 
fT'„  Thermal  residual  stress  in  the  whisker 
Au  Difference  in  the  coefficients  of  thermal  expansion 
A7  Magnitude  of  temperature  change 
Ah  Difference  of  displacement  at  the  w  hisker  matrix 
interface  and  longitudinal  displacement  away  from 
the  interface 
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t  lnlerlaci.il  shear  stress 
Ait.  An  when  matrix  apptoaches  welding  la  ,  ) 

/)  and  d  Diameters  ot  the  unit  cell  and  whisker 
l„  Whisker  volume  traction 

I ...  Young  s  modulus  ol  the  matrix 
/■..,  Y  oung's  modulus  ol  the  whtskets 

1.. ..,  Work-hardening  rate  of  the  matrix 
o  ...  Matrix  \  icld  stress 

/  Total  length  of  the  whisker 

Iaiigeiiti.il  shear  modulus  ot  the  matrix 

tr  Applied  stress 

1 .. ..  Total  applied  strain 

Average  longitudinal  matrix  lesidual  stress  ,tv  ., 
result  ol  external  applied  stiess  ,r 
I.  .  Young's  modulus  ot  the  composites 
(  T.  and  (  Stillness  and  tangential  stitfne"  of  the  matrix 
r  and  it  Average  sti.un  and  stress  on  the  weld  surfaces  ill 
the  strain  and  stress  space 
i',  Poisson's  ratio  of  the  constituent  in  the  composites 
< •  r ‘  Average  longitudinal  thermal  iesidu.il  stress  m 
the  matrix 

V '...,  Volume  fraction  ot  tile  initial  tension  /one  in  the 
matrix 

or;.,.  Average  longitudinal  matrix  thermal  residual  stress 
in  the  initial  tensile  /one 

ir...  i  Average  thermal  residual  stress  in  the  compics- 
sivc  /one 

A /  Magnitude  of  the  temperature  change  that  will 
induce  plastic  deformation  at  the  top  ot  the 
vvhiskci 

A 7,  Magnitude  of  tile  temperature  change  b\  which 
the  entire  matrix  is  plastically  deformed 
'.  '<K,  >  Average  long.tudin.il  internal  stress  m  the  initial 
tensile  and  compressive  /one  respectively .  when 
a  uniaxial  external  applied  stress  is  applied 
ur (<r;.,;i)  Average  longitudinal  residual  stress  m  the  initial 
tensile  and  compressive  /one  after  the  uniaxial 
applied  external  load  is  released 
t  Applied  plastic  strain 

Introduction 

A  previous  investigation  of  the  deformation  behavior  of  short 
whisker-reinforced  silicon  carbide-aluminum  (Nit '-All  compos 
lies  |/|  has  indicated  that  the  constitute  behavior  of  these  com¬ 
posites  is  distinctly  different  under  uniaxial  compressive  defor¬ 
mation  as  compared  to  tensile  deformation  It  was  generally 
found  that  the  compressive  weld  strength  of  these  composites 
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was  larger  than  that  of  the  tensile  yield  strength,  while  the  ap¬ 
parent  Young's  modulus  for  compression  was  smaller  than  that 
in  tension.  If  monolithic  aluminum  is  initially  deformed  plasti¬ 
cally  in  tension,  then  reverse  deformed  in  compression,  a  dif¬ 
ference  in  stress  is  required  to  initiate  plastic  deformation  in  the 
reverse  cycle.  Usually,  the  magnitude  of  the  yield  stress  in  the 
reverse  cycle  is  less  than  the  flow  stress  in  the  forward  cycle. 
This  reduction  in  stress  (or  the  strain  necessary  to  reach  the 
previous  level  of  flow  stress)  is  defined  as  the  Bauschinger  effect. 
In  the  case  of  monolithic  aluminum,  the  direction  of  initial  de¬ 
formation  (that  is.  tensile  or  compressive)  does  not  have  any 
influence  on  the  Bauschinger  effect.  The  SiC-AI  composites, 
however,  exhibit  a  remarkable  characteristic,  that  is.  the  flow 
stress  drop  (or  alternatively,  the  Bauschinger  strain)  is  larger 
when  the  composite  is  first  precompressed  plastically  and  sub¬ 
sequently  deformed  in  tension  than  when  the  composite  is  first 
deformed  in  tension  and  then  in  compression  |/J. 

It  is  the  opinion  of  the  authors  that  the  aforementioned  dif¬ 
ferences  in  the  constitutive  behavior  under  tensile  and  com¬ 
pressive  loading  as  well  as  the  asymmetric  Bauschinger  phenom¬ 
enon  are  primarily  due  to  the  influence  of  the  residual  stresses. 
It  has  been  shown  that  thermal  residual  stresses  affect  not  only 
the  elasto-plastic  deformation  behavior  of  the  composites  [ /  ] . 
but  also  the  fracture  behavior  [2|.  Therefore,  it  becomes  increas¬ 
ingly  important  to  investigate  the  influence  of  residual  stress  on 
the  deformation  process 

The  material  to  be  simulated  in  this  investigation  consists  of 
an  annealed  6061  aluminum  alloy  matrix  reinforced  by  silicon 
carbide  whiskers.  The  primary  goal  of  this  paper  is  to.  first, 
investigate  analytically  the  generation,  arrangement,  and  mag¬ 
nitude  of  the  residual  stresses  caused  by  the  difference  in  coef¬ 
ficients  of  thermal  expansion  (ACTE)  between  the  reinforcement 
and  the  matrix  during  the  cooling  process;  second,  studs  the 
changes  in  the  residual  stress  as  a  result  of  external  deformation 
either  in  tension  or  compression;  and  finally,  investigate  the 
differences  in  the  constitutive  behavior  under  different  external 
loading,  that  is.  deformation  conditions. 

In  an  experimental  procedure  described  elsewhere  ( /..•?(.  the 
materials  is  annealed  at  723K  and  then  tested  either  in  uniaxial 
tension  or  compression  at  room  temperature.  Our  commitment 
in  this  investigation  is  to  adhere  to  this  physical  process  and. 
thereby,  to  simulate  the  constitutive  behavior  of  the  composites 
and  monitor  the  change  of  the  residual  stresses  These  residual 
stresses  were  induced  during  the  cooling  process  and  subse¬ 
quently  altered  by  external  load 

While  several  investigators  have  concentrated  their  efforts  on 
the  thermal  residual  stresses  [J-<Y|.  others  [V  !S\  have  worked 
on  the  elasto-plastic  behavior  of  the  short  whisker-reinforced 
composites  w  ith  emphasis  on  the  stress-strain  relationship  of  the 
composites,  that  is.  the  constitutive  behavior.  Their  research 
includes  numerical  treatments  that  use  finite  element  procedures 
[9-/.?].  Some  of  them  paid  special  attention  to  the  microstruc- 
tural  change  in  the  matrix  \I0.I. I |.  while  others  were  based  on 
various  theoretical  bases  with  different  assumptions  |/4-  IS |. 

Recently.  Arsenault  and  Taya  |/9|  and  Withers  et  al  (2(t|  used 
the  Tshelby  method  |2/.22|  to  explain  the  yielding  behavior  of 
the  short  whisker-reinforced  metal-matrix  composites  (MMC) 
under  the  influence  of  the  thermal  residual  stresses.  They  found 
that  the  residual  stresses  are  totally  responsible  for  the  difference 
in  the  yield  strength  of  the  whisker-reinforced  MMC  between 
tensile  and  compressive  uniaxial  loading.  However,  despite  the 


effort  by  the  previous  investigators  [19,20 J.  the  role  residual 
stresses  played  in  the  constitutive  behavior  of  the  composites  or 
the  effect  of  the  external  loading  on  the  formation  of  the  residual 
stresses  has  never  been  systematically  investigated.  Nevertheless, 
these  experimental  and  theoretical  results  1 79,20)  have  provided 
us  with  a  clear  indication  that  the  residual  stresses  may  be  a 
dominating  factor  that  accounts  for  the  differences  of  the  con¬ 
stitutive  behavior  of  the  composites  under  a  uniaxial  tensile  or 
compressive  load. 

The  analytical  investigation  in  this  study  is  carried  out  by  two 
distinct  methods:  a  finite  element  analysis  and  a  phenomeno¬ 
logical  modelling  scheme  involving  two  models.  The  main  ob¬ 
jective  of  the  finite  element  analysis  is  to  monitor  the  stress- 
strain  relationship  and  the  changes  that  occur  in  the  value  of  the 
residual  stresses  as  a  result  of  the  external  deformation.  The 
purpose  of  the  phenomenological  models,  on  the  other  hand,  is 
to  obtain  some  physical  insight  into  the  various  factors  that  con¬ 
trol  the  processes  of  residual  stress  generation  and  alteration  and 
the  role  these  play  in  influencing  the  constitutive  behavior  of  the 
composites. 


The  Finite  Element  Model 

To  determine  the  response  of  the  composite  material  under 
the  influence  of  residual  stresses,  a  simplified  two-dimensional 
(2-D)  plane  strain  finite-elem  analysis  was  carried  out  using 
the  ADINA  code.  In  this  analysis,  the  reinforcement  is  assumed 
to  be  perfectly  elastic,  the  matrix  is  assumed  to  follow  bilinear 
stress-strain  relationship,  and  the  isotropic  hardening  rule  is  used 
Particular  interest  was  paid  to  determining  the  macroscopic  stress- 
strain  behavior  of  the  composite  and  performing  a  microscopic 
analvsis  of  the  subsequent  changes  that  occur  in  the  values  of 
the  matrix  longitudinal  residual  stresses  after  the  external  load 
was  removed.  The  modelling  of  the  change  'f  the  residual  stress 
as  a  function  of  applied  load  was  also  performed  through  three- 
dimensional  (3-D)  finite  element  method  (FFM)  analysis  with 
the  same  boundary  condition  consideration.  The  results  are  com¬ 
pared  with  those  obtained  from  the  2-D  plane  strain  model. 

A  recent  studv  h\  Takao  and  Taya  [3T|  showed  that  the  use 
of  the  volume-average  aspect  ratio  will  provide  an  accurate  ac¬ 
count  for  the  overall  properties  ol  the  composites,  and  the  em¬ 
ployment  of  the  perfect  alignment  assumption  will  give  rise  to 
satisfactory  results  Therefore,  in  the  analyses,  the  value  of  the 
volume  fraction  of  the  whiskers  as  well  as  the  whisker  aspect 
ratio  were  held  constant  (V„  26'r  and  Id  4)  A  value  of 

Id  4  is  a  typical  representation  of  the  average  value  observed 
experimentally  |24.J.,'|  Furthermore,  it  was  assumed  that  the 
composite  consists  of  an  array  of  hexagon. illy  distributed  par¬ 
allelepiped  SiC '  whiskers  embedded  m  a  6061  aluminum  alloy 
matrix  (Fig  I )  A  unit  cell  can  be  constructed  because  of  periodic 
reinforcement  arrangement.  Previous  investigations  |9.//|  and 
our  current  analysis  have  indicated  that  within  a  particular  re¬ 
gion.  a  small  variation  in  the  values  of  interparticle  spacing  does 
not  significantly  affect  the  material  behavior  (tor  example,  tor 
vield  stress,  see  Fig.  2).  The  value  of  l  H  in  Fig  1.  therefore, 
was  not  selected  as  a  critical  parameter  in  the  analysis  and  was 
kept  as  a  constant  (11.29)  throughout  the  analyses  Note  that 
although  from  Fig  1  it  appears  that  the  effect  of  the  variation 
of  V  H  is  slightly  smaller  when  a  larger  V  II  value  is  selected, 
the  rate  of  convergence  is  better  when  the  current  value  was 
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FIG  1  —  The  whiskers  are  distributed  as  a  periodic  liesaitonat  infinite 
array  in  the  matrix  flv  considering  symmetry,  a  unit  cell  that  contains 
only  a  quadrant  oj  the  particle  is  needed,  as  indicated  in  the  future 


0  0  3.0  6  C 

V/H 

FIG.  1  —  HIM  read:,  nniuult  dial  llic  ' itid  stress  ol  tin  i nmposucs  is 
not  sensitis  e  to  the  ratio  of  V  II  defined  in  lie  I 

selected,  that  is.  it  will  result  in  a  more  efficient  usage  of  com¬ 
puter  time. 

Figure  3  indicates  the  exact  details  of  the  finite-element  mesh 
Only  a  quadrant  of  the  composite  material  which  contains  the 
SiC  whisker  was  selected  to  represent  the  whole  infinite  array 
of  fibers  in  the  matrix  A  multi-constraint  boundary  condition 
was  used  in  the  analysis  to  consider  the  geometrical  symmetry 
and  the  condition  of  the  internal  stress  indicated  as  the  following: 

|  <r , ,n, til  ~  <t' 

where  <r„  is  the  stresses  at  the  boundary.  ir„'  is  the  applied  stress. 
n,  is  the  unit  outward  normal,  and  /,  represents  four  different 
sides  of  the  boundary  of  the  unit  cell. 

Since  the  constitutive  behavior  of  the  SiC  Al  whisker  com¬ 
posites  has  been  extensively  investigated  experimentally  |/.35.3ft|. 


FIG  hi  \1  mesh  is  here  '/«•  to  entht  noile  elements  ssere  used 


it  is  appropriate  to  compare  our  FF.M  results  with  the  experi¬ 
mental  results  The  experimental  evidence  |/|  indicated  that  the 
apparent  Young's  modulus  and  the  yield  strength  ot  the  com¬ 
posites  were  different  when  a  uniaxial  tensile  and  compressive 
load  was  applied.  The  apparent  Young's  modulus  is  larger  under 
tensile  uniaxial  applied  load,  whereas  the  yield  strength  is  higher 
under  compressive  load  It  is  a  hypothesis  at  this  stage  that  the 
differences  in  the  constitutive  behavior  ot  the  composites  are  due 
mainly  to  the  presence  of  the  thermal  residual  stresses  There¬ 
fore.  by  considering  the  initial  thermal  residual  stresses  which 
are  generated  during  the  cooling  process,  such  differences  should 
be  reproduced  bv  the  FI  M  simulation  Consequently,  it  such 
differences  are  picdicted.  then  the  residual  stiesses  that  are  te- 
sponsible  tor  producing  such  differences  should  be  reasonably 
well  predicted.  Figure  4  shows  the  uniaxial  stress-strain  curves 
produced  by  the  FF.M  analysts  where  the  global  stress  and  strain 
are  defined  as: 

it 

t  - :  it  or ) 

T  -  / 

where  it,  and  (o,>  are  the  average  displacement  and  normal  stress 
along  the  boundary  whose  normal  is  parallel  to  the  direction  ot 
the  uniaxial  loading.  The  value  mi  can  also  be  defined  as  the 
applied  load  intensitv  (</)  and  \  .l  are  defined  in  I  ig  1 

Comparing  the  stress-strain  relationship  under  tension  and 
compression,  it  can  be  seen  that  the  differences  in  Young's  mod¬ 
ulus  and  yield  strength  of  the  composites  have  been  reproduced 
in  a  correct  wav 

A  Hauschingcr  test  can  also  be  carried  out  bv  I  I  M  I  he  di¬ 
rection  of  initial  loading  has  an  effect  on  the  magnitude  ol  the 
Hauschingcr  strain  (i, ).  as  defined  ill  Fig  3  It  initiallv  tested  in 
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FIG.  4— Stress-strain  curve s  predicted  hy  FEM  » here  die  apparent 
Young  v  modulus  of  the  composite  is  higher  a  lien  it  is  in  tensile  toad, 
whereas  the  yield  stress  is  greater  a  hen  the  sample  is  under  compressive 
load 


FIG  5  —  FEM  prediction  of  Bauschtnger  strain  h  h>  where  t is  larger 
if  the  sample  is  under  the  compression- first  loading  scheme 

tension  and  then  reverse  loaded  into  compression,  the  Baus- 
chinger  strain  is  smaller  compared  to  the  case  in  which  initial 
loading  is  in  compression  followed  hy  tension  (Fig.  5).  This  has 
been  observed  experimentally  [/|. 

In  the  study  of  the  internal  stresses  in  composite  materials,  it 
is  a  common  practice  to  use  the  concept  of  volume  average  both 
experimentally  (J7. 2<S]  and  theoretically  (/9, 20.29],  In  this  re¬ 
search.  the  same  method  is  employed.  Because  the  longitudinal 
thermal  residual  stress  is  much  larger  than  the  transverse  residua! 
stress  [/9 j.  it  is  believed  that  the  observed  differences  in  consti¬ 
tutive  behavior  are  due  mainly  to  the  longitudinal  thermal  re¬ 
sidual  stress.  The  longitudinal  thermal  residual  stress  in  the  com¬ 
posite  was  therefore  determined  by  first  evaluating  the  residual 
stress  in  finite-element  mesh  and  then  averaging  these  values 
over  the  whole  specimen.  The  overall  average  residual  stress  in 
the  matrix  was  calculated  through  a  simple  averaging  scheme 


given  by  the  following  equation: 

£  KhA 

<«V  =  ‘  y  —  (/  =  j  =  3)  (3) 

k 

where  (cr,,)4  is  the  stress  in  element  k ,  and  /14  is  the  area  of  that 
element. 

Using  such  a  simplified  averaging  scheme,  it  was  found  that 
upon  application  of  1  c'r  plastic  strain,  if  the  applied  loading  is 
tensile,  the  average  matrix  residual  stress  changes  from  its  initial 
tensile  state  of  47.5  MPa  to  the  compressive  state  of  -  10.8  MPa. 
while  if  the  applied  loading  is  compressive,  the  average  matrix 
longitudinal  (tensile)  residual  stress  drops  20.1  MPa  and  then 
increases  back  up  to  29.3  MPa  when  it  reache-s  1  rr  plastic  strain 
The  overall  change  in  residual  stress  as  a  function  of  applied 
strain  is  shown  in  Fig.  6.  It  is  clear  that  the  primary  effect  on 
the  longitudinal  residual  stress  from  compressive  loading  is  that 
the  magnitude  of  the  residual  stress  in  the  matrix  reduces  and 
then  increases  upon  further  loading  while  tensile  deformation 
produces  an  overall  reduction  in  the  tensile  residual  stress  and 
eventually  changes  the  residual  stress  from  tensile  to  compres¬ 
sive.  That  is.  the  change  in  residual  stress  does  not  occur  mon- 
otonically  with  increasing  deformation. 

ft  is  also  interesting  to  see  that,  in  the  case  of  compressive 
deformation,  the  initial  tensile  residual  stress  reduces  until  the 
total  applied  strain  reaches  about  0.3G  which  is  within  close 
proximity  of  0.2 rr  strain  of  yielding.  At  this  point,  the  tensile 
residual  stress  reaches  a  minimum.  Upon  further  deformation, 
there  is  an  increase  in  the  tensile  residual  stress.  In  the  case  of 
tensile  deformation,  the  tensile  residual  stress  initially  decreases 
at  a  more  rapid  rate  to  -0.3G  total  strain,  then  the  rate  of 
decrease  becomes  less.  At  a  total  strain  of  — O.KG .  the  residual 
stress  becomes  compressive.  The  longitudinal  residual  stress  con¬ 
tour  after  different  levels  of  the  total  applied  strain  is  shown  in 
Figs.  7  through  1 1 . 

Despite  the  fact  that  the  total  average  longitudinal  thermal 
residual  stress  is  in  tension,  previous  FEM  results  [/.?]  have  in- 


FIG  ft— .4  verage  longitudinal  residual  stress  as  a  fum  turn  of  the  total 
applied  strain  preda  ted  by  bF.M 
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FIG.  1 1  — Longitudinal  residual  stress  contour  at  a  total  applied  strain 
of  =  IM'i 


dicated  that  in  whisker-reinforced  composite  materials,  the  re¬ 
sidual  stress  distribution  is  such  that  near  the  tip  of  the  whisker, 
the  matrix  is  under  compression  while  the  rest  of  the  matrix  is 
in  tension,  thus  giving  rise  to  an  overall  tensile  average  residual 
stress  in  the  matrix.  Interestingly  enough,  our  finite  element 
analysis  indicated  that  under  applied  tensile  load,  the  compres¬ 
sive  and  tensile  zones  exchanged  their  signs  within  a  reasonable 
range  of  strain  (Figs.  7-11).  while  compressive  loading  led  to  a 
residual  stress  increase  in  both  regions.  Figure  12  is  a  schematic 
diagram  of  the  tensile  deformation-induced  redistribution  of  the 
longitudinal  thermal  residual  stress. 

Another  persuasive  fact  to  be  noted  is  that  while  the  global 
average  residual  stress  exhibited  a  non-monotonic  behavior,  the 
residual  stress  near  the  whisker  matrix  interface  never  showed 
an  initial  drop  after  loading  as  it  did  in  the  case  of  global  average 
residual  stress.  The  response  of  the  residual  stress  to  the  external 
load  in  the  interface  was  such  that  the  longitudinal  residual  stress 
increased  upon  an  initial  compressive  plastic  prestrain  and  de¬ 
creased  with  an  initial  tensile  plastic  prestrain,  as  can  be  seen  in 
Fig.  13. 

Although  the  simplified  averaging  scheme  discussed  above  is 
easy  to  use.  it  reflects  the  nature  of  the  stress  distribution  of  a 
uniform  field.  It  is  evident  that  this  scheme  does  not  reflect  the 
fact  that  the  deformation  upon  loading  is  inhomogeneous  in  the 
matrix.  Because  of  this  deformation  inhomogeneity,  for  the  same 
applied  load  increment,  the  importance  of  the  residual  stress  will 
vary  depending  on  the  local  stress  concentration.  The  higher  the 
load  deformation  rate,  the  higher  the  stress  concentration.  Fur¬ 
thermore.  the  difference  in  the  value  of  the  local  defoimation 
rate  at  different  locations  in  the  matrix  constitutes  the  defor¬ 
mation  gradient  in  the  matrix.  It  is  well  known  that  a  gradient 
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FIG.  12 — As  a  result  of  tensile  loading,  the  arrangement  of  residual 
stress  alters  so  that,  at  the  whisker  up.  the  sign  of  the  residual  stress  changes 
from  tension  to  compression. 
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FIG.  1.1 — Average  residual  stress  near  the  longitudinal  matrix  whisker 
interface  as  a  function  of  applied  stress  The  residual  stress  is  taken  from 
an  area  of  Od I  *  0.2 d  around  the  transverse  centerline  of  the  whisker 
near  the  interface 


in  plastic  deformation  is  responsible  for  the  deformation  gen¬ 
erated  residual  stress. 

In  the  previous  investigations,  it  was  a  common  practice  to 
assume  that  residual  stress  was  generated  due  only  to  the  de¬ 
formation  gradient  between  the  reinforcement  and  the  matrix 
\20\.  Therefore,  to  evaluate  best  the  effect  of  the  thermal  residual 
stress  on  the  subsequent  loading  and  the  degree  of  matrix  de¬ 
formation  inhomogeneity,  a  new  averaging  rule  should  be  de¬ 
veloped.  The  authors  believe  that  application  of  a  weighted- 
averaging  scheme  is  more  appropriate  and  will  lead  to  a  more 
accurate  account  of  the  actual  physical  process.  That  is.  to  in- 
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corporate  the  local  loading  rate  (dcrJUq),  as  a  measure  of  defor¬ 
mation  gradient,  directly  into  the  average  stress  formula. 

An  averaging  formula  of  the  type 


where 

(Acr„)  -  -  (c7„)J,  (4 b) 

and 

■*</  =  </l, .  a,  -  </l,  (4<  ) 

was  then  used,  where  q  is  the  applied  load  intensity  and  <7,(  now 
denotes  the  new  average  residual  stress.  From  Eq  4a  it  can  be 
seen  that  for  homogeneous  matrix  deformation,  the  stress  rate 
(Acr./Aq)*  is  a  constant  over  the  entire  matrix.  Therefore,  the 
average  residual  stress  predicted  by  Eq  4 a  is  the  same  as  the 
prediction  made  through  Eq  3.  Evaluating  the  average  matrix 
residual  stress  using  Eq  4rr.  it  was  found  that  upon  an  initial  load, 
the  ratio  between  the  two  averages  is: 

^  -  1.35  (i  -  I  -  3)  (5) 

,T„ 

The  factor  1.55  is  an  indication  of  the  degree  of  deformation 
heterogeneity  of  the  composite  material.  liquation  5  indicates 
that  the  simplified  average  thermal  residual  stress  is  greater  than 
the  one  produced  from  the  weighted  scheme,  which  means  that 
the  weighting  factor  is  larger  in  the  region  where  the  residual 
stress  has  a  negative  value.  This  simply  means  that  deformation 
in  the  compressive  zone  (tip  of  whiskers)  is  more  severe  with  a 
liaxial  applied  external  load. 


C 


B 

A 

IK  i  14  II.  \1  results  militated  that  the  residual  strew  field  tan  he 
characterized  mio  three  regions  m  terms  nf  strew  rate,  w  here  at  the  w  luskei 
ttp.  the  stress  rate  is  highest 


To  understand  better  the  heterogeneous  deformation  field  in 
the  matrix  as  a  result  of  the  existence  of  the  w  hiskers,  an  average 
stress  rate  (da.liq)  diagram  was  constructed  for  the  composite  at 
different  points  of  the  matrix  when  the  composite  is  initially 
loaded  in  tension,  as  shown  in  Fig.  14.  In  particular,  three  regions 
within  the  matrix  were  of  interest.  The  analysis  indicated  that 
the  average  stress  rate  within  the  Regions  A.  B.  and  C  were, 
respectively.  6  x  10  -.9  x  Id  -'.and  13  x  111  (MPa  MPa) 
These  results  support  the  notion  that  the  compressive  zone  at 
the  tip  of  the  whisker  responds  faster  to  external  loading  than 
any  other  region  in  the  matrix.  On  the  other  hand,  in  the  vicinity 
of  the  longitudinal  whisker-matrix  interface,  the  response  of  the 
internal  stress  to  applied  loading  is  the  slowest.  Finally,  in  Region 
b.  where  some  interaction  with  the  whisker  tip  is  expected,  the 
response  to  applied  loading  falls  within  an  intermediate  state. 

The  finding  that  there  is  a  stress  concentration  at  the  whiskers 
tips  is  consistent  with  what  was  observed  experimentally.  Sato 
et  al.  [30]  have  reported  extensive  damages  at  short  fiber  tips  as 
a  result  of  excess  plastic  deformation  Note  that  the  low  defor 
mation  rate  in  Region  A  is  due  mainly  to  the  interfacial  load 
transfer,  whereas,  the  disparity  of  deformation  rate  between  Re¬ 
gions  B  and  C  is  a  result  of  be  volume  mismatch  between  the 
reinforcement  and  the  matr  > 

In  consideration  of  the  deformation  process  in  the  matrix,  it 
is  believed  that  the  differences  of  microstrain  rate  in  different 
parts  of  the  matrix,  that  is.  inhomogeneous  matrix  deformation, 
are  solely  due  to  the  presence  of  the  whisker.  That  is.  near  the 
longitudinal  whisker-matrix  interface,  deformation  is  restricted 
because  of  the  high  stiffness  of  the  SiC  whisker,  and  at  the  same 
time,  the  matrix  material  at  the  whisker  tip  deforms  more  to 
accommodate  the  global  strain.  Therefore,  the  nature  of  the 
deformation  gradient  is  not  expected  to  exhibit  a  large  change 
during  the  entire  course  of  deformation 

The  finite-element  analysis  presented  in  this  paper  has  been 
performed  with  the  sole  intention  of  obtaining  a  better  under¬ 
standing  of  the  effect  of  the  residual  stresses  on  the  constitutive 
behavior  of  discontinuouslv  reinforced  metal-matrix  composites. 
It  is  evident  that  the  most  prominent  characteristics  of  the  con¬ 
stitutive  behavior  of  the  composites  have  been  adequately  pre¬ 
dicted  using  2-D  FEM.  The  physical  reason  tor  the  effectiveness 
of  the  current  2-1)  analysis  is.  as  shall  be  seen  in  the  later  dis¬ 
cussion.  that  the  2-1)  composite  does  not  destroy  the  inherent 
mechanism  of  residual  stress  generation  However,  note  that 
although  the  2-1)  FEM  model  provides  a  good  economical  tool 
for  understanding  how  the  system  operates,  it  was  not  the  in¬ 
tention  ol  the  authors  in  the  present  investigation  to  compare 
quantitatively  the  analytical  results  with  those  obtained  experi- 
mentallv  It  is  evident  that  the  plane-strain  representation  ol  an 
otherwise  much  more  complex  3-1)  array  ot  whiskers  in  a  matrix 
is  an  oversimplification  toi  quantitative  analysis.  Figure  15  shows 
a  comparison  of  residual  stress  as  a  function  ol  applied  strain 
predicted  by  the  plane  strain  model  and  the  actual  31)  FEM 
with  the  same  volume  traction  Note  that  there  is  no  fundamental 
difference  in  the  pattern  ot  the  change  in  residual  stress  as  a 
result  of  external  deformation  The  only  feature  noticeable  is  the 
steeper  overall  slope  predicted  by  the  3-1)  I  I  M 

Phenomenological  Modelling 

Ihe  results  ot  the  finite-element  analvsis  discussed  in  the  pre¬ 
vious  section  indicate  that  the  magnitude  ot  the  residual  stress 
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FIG.  15  —  Comparison  of  average  residual  strew  preda  ted  b\  2-1)  and 
2-1)  H.M 

as  a  result  of  ACTE  and  its  subsequent  change  under  the  influ¬ 
ence  of  mechanical  deformation  are  primarily  governed  by  the 
deformation  gradient  in  the  matrix.  It  is  our  belief  that  the  ex¬ 
istence  of  the  initial  compressive  zone  is  a  direct  consequence  of 
volume  mismatch  between  the  whisker  and  the  matrix,  which 
occurs  as  a  result  of  ACTE  in  the  composite.  On  the  other  hand, 
the  fact  that  the  average  matrix  residual  stress  is  in  tension  can 
also  be  obtained  through  interface  shear  load  transfer.  It  is. 
therefore,  clear  that  the  process  of  residual  stress  generation  and 
residual  stress  change  are  basically  controlled  by  two  factors:  the 
nature  of  the  interfacial  shear  load  transfer  between  the  SiC 
whiskers  and  the  Al  matrix  and  the  volume  mismatch  between 
the  two  constituents.  The  mismatch  of  the  volume  is  introduced 
by  either  the  CTE  mismatch  when  under  thermal  loading  or  the 
differences  in  Young's  modulus  and  Poisson's  ratio  when  under 
external  loading. 

To  maintain  overall  compatibility,  the  two  factors  must  work 
together.  The  volume  difference  is  matched  through  interfacial 
normal  stress,  and  the  interfacial  compatibility  is  satisfied  by 
means  of  interfacial  shear  stress.  To  comprehend  the  separate 
roles  of  "interface  effect"  and  "volume  mismatch"  on  the  resid¬ 
ual  stress  response  after  the  mechanical  loading,  two  distinct 
analytical  models  were  developed.  The  first  model,  subsequently- 
referred  to  as  the  "interface  model."  is  shown  in  Fig.  16,  which 
assumes  that  the  material  consists  of  a  discontinuous  filament 
embedded  in  the  matrix.  It  is  essentially  the  same  as  the  shear 
lag  model  with  a  consideration  of  thermal  stress  and  average- 
plastic  deformation.  It  is  also  assumed  that  there  is  full  elastic- 
continuity  across  the  whisker-matrix  interface.  The  whiskers, 
therefore,  remain  fully  bonded  to  the  matrix,  --ausing  complete- 
load  transfer  between  the  two  constituents.  The  primary  purpose 
of  this  type  of  modelling  is  to  understand  the  effect  of  load 
transfer  at  the  interface  on  the  residual  stress  as  a  result  of  the 
applied  loading. 

In  the  second  model,  it  is  assumed  that  a  discontinuous  whisker 
is  embedded  in  the  matrix.  To  discount  the  effect  of  interfacial 
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FlCi  16 — Schematic  diagram  lor  the  interline  model 

shear  load  transfer,  interfacia!  sliding  is  allowed  across  the  whis¬ 
ker-matrix  interface  along  the  direction  of  the  whisker  axis,  but 
there  is  full  continuity  of  the  interfacial  normal  stress  and  dis¬ 
placement  at  the  tip  of  the-  whisker  between  the  whisker  and  the 
surrounding  matrix.  It  is  further  assumed  that  a  boundary  is  made 
along  the  extension  of  the  longitudinal  Sit'  Al  interface,  as  shown 
in  Fig.  17.  This  model,  which  is  referred  to  as  the  "volume- 
mismatch  model."  is  show  n  in  Fig.  17.  The  main  purpose  of  this 
model  is  to  explain  how  the  volume  mismatch  affects  the  residual 
stress  arrangement  as  a  result  of  a  temperature  change  and  ex¬ 
ternal  loading. 

In  both  of  these-  analyses,  the  matrix  material  is  assumed  to 
follow  a  bilinear  stress-strain  relationship,  and  the  reinforcement 
is  assumed  to  be  perfectly  elastic.  The  two  analy  tical  models  are- 
discussed  in  more  detail  below. 

The  Shear  Lag  Model 

The  simplest  and  yet  most  effective  approach  to  account  lor 
the  mechanism  of  shear  load  transfer  at  the  matrix-reinforcement 
interface  is  the*  "shear  lag"  theory  introduced  by  Cox  [7/|.  Ever 
since  the  advent  of  this  theory  ,  it  has  been  widely  used  because 
of  its  simplicity.  Some  investigators  [.?_?. 77|  have  claimed  ii  to  be 
quite  accurate  in  predicting  mechanical  properties  of  composites; 
others  |25..T/|  have  found  it  less  attractive.  In  this  section  a 
modified  "shear  lag"  method  is  used  to  try  to  account  for  the 
generation  of  the  residual  stresses.  It  is  assumed  that  the  overall 
constitutive  response  of  the  matrix  remains  elastic  until  the  av¬ 
erage  matrix  stress  state  becomes  plastic. 

When  the  composite  is  cooled  down  from  the  annealing  tem¬ 
perature.  longitudinal  thermal  residual  stresses  are  developed  in 
the  matrix  and  the  whisker.  These  will  be  denoted  as 
(irL.)  and  <<r'„).  respectively.  Using  a  modified  shear  lag  model 
that  incorporates  temperature  effect  and  plasticity  ,  appropriate- 
expressions  may  be  readily  derived  for  (<0  and  The  gov¬ 
erning  equation  for  the  modified  shear  lag  model  is  (further 
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FICi.  17  --  Vi  ItcHiiitn  i havr<mi  lor  the  voliwn  niivniilt  h  model.  where 
inlcrtueuil  bonding  n  iiinored 


stress.  /  is  the  total  length  ot  the  whisker,  .tnd  h  and  (i  can  he 
determined  as  follows: 


NO',,,., 

d  Itlt  I)  d) 


-  \  ••  (>,.,,1  . 
d  \  ln(  l)  1 1 ) 


IN  />) 


where  O',,,.,  is  the  tangential  shear  modulus  of  the  matrix 

When  a  uniaxial  loading  hr  )  corresponding  to  the  total  strain 
(r,.„ )  is  applied  along  the  whisker's  direction  in  either  tension  or 
compression,  and  the  composite  is  then  unloaded,  it  max  be 
shown  that  the  matrix  residual  stress  (ur:.;>)  as  a  result  ot  applied 
loading  now  becomes: 


l A. . ) 

t;,"  '  -  or;,.i  - r 

A.  , 


(  Hint 


details  are  given  in  Appendix  A  and  elsewhere  |.’’5|): 

dir,,  2D  4t 

——  /it  An!  A  /  c  -  Aic  -  An  )  —  ir  -  —  (01 

d:  d  d 

where  A<»  is  ACTF;  A  /  is  the  magnitude  ot  temperature  change. 
Aic  -  n  -  t  where  u  and  i  arc  the  displacements  along  whiskei 
axes  at  and  away  from  the  interface  in  the  matrix.  An  represents 
Aw  when  the  matrix  becomes  plastically  deformed.  <r  „  is  the 
matrix  yield  stress.  ir„  is  the  whisker  stress,  li  is  a  constant  to  be 
determined,  and  I)  and  d  are  the  diameters  of  the  unit  cell  and 
the  whisker,  respectively.  It  the  temperature  differential  (  A/  ) 
is  such  that  the  strains  in  both  the  matrix  and  the  whisker  lemam 
elastic,  then  Aw  and  ir  ..  may  be  taken  away  Iron)  I  q  o.  and 
Fq  b  will  reduce  to  the  form  given  In  Cox  |.*/|  with  a  new  thermal 
mismatch  term  t/i  A«r  A  /  )  For  a  2()A  Ni(  Al  composite,  it  is 
likely  that  the  matrix  will  be  plastically  deformed  as  a  result  ot 
the  cooling  process.  Following  the  procedure  summarized  by 
Kelts  |./b|.  the  general  form  of  the  solution  is  obtained  as  follows 
(Sec  Appendix  A): 


tir;„i 

I  *  r.  tanli  ( [il  2) 

*  I  V  .  A  ..."  (|i/2) 

(Ait 

or,,  > V  „  •  <ir;„  >(  I  l  , )  II  1 7/i| 


where  V’„  is  the  volume  fraction  ot  the  whiskers.  and  A... 
represent  the  Young's  moduli  of  the  matrix  and  whisker.  A..,,,  is 
the  work-hardening  rate  of  the  matrix.  <r,„,  is  the  matrix  yield 


tanli  I  ft/  2 1 

nn  i)  ~ 


Ac.  A/ 


l  . 

or. -i  -  -  •  <r;:;  •  till/') 

:\ :  /. 


where  i  .,  represents  the  applied  plastic  strain  and  /.  is  the 
Young's  modulus  of  the  composites.  Here,  as  m  the  oriental 
references  xl  and  the  1  Cung's  modulus  ot  the  composites  is 
assumed  to  tollow  the  rule  ot  mixture  | A.  A.,,(  1  f  .i  - 
A...  I  „|  Smcc  A.,,,,  clearly .  the  matrix  residual  stress  de¬ 

creases  upon  tensile  loading  tr  ,,  III  and  increases  on  com¬ 
pressive  loading  it  ,.  0)  Flic  overall  changes  in  the  residual 

stress  remain  monotonic 

I  lie  sel  ot  equations  state  above  appiv  equallv  to  both  cases 
ot  compicssive  and  tensile  loading  provided  we  considei  '.lie 
trailslormation: 


1  his  means  that  cry  is  svmmctric  about  cr;  when  plastic 
dctormalion  is  coiisideicd  It  also  can  be  seen  that  Fq  Id  is  only 
a  function  ot  applied  plastic  strain,  vv  Inch  means  that  clastic  strain 
will  not  alter  the  value  ot  the  residual  stress  which  contradicts 
the  fact  that  there  is  ,m  initial  icsidual  stiess  drop  regardless  of 
the  loading  directions  predicted  In  FI  \l  It  is.  thcrctoic.  deal 
that  the  modified  shear  lag  model  cannot  adequately  pi  edict  the 
effect  ot  the  loadin'!  directions  on  the  usvmmctric  behavioi  ot 
the  average  residual  stiess  It  otliet  wools,  the  mechanism  ot 
load  transfer  throuiih  the  matrix-w hiskei  interlace  cannot  have 
a  predominant  effect  on  the  aveiage  value  ot  the  icsidual  stress 
and  the  residual  strese  chance  and.  coitscqucntlv .  on  the  con 
st  it  li  t  iv  e  behavioi  In  tact,  the  outcome  ot  tills  model  matches 
the  results  ot  the  tunic  element  unulvsis  eivcn  earlier  that  within 
the  immediate  proxmutv  ot  the  longitudinal  mutiix  reinforce¬ 
ment  interface,  the  icsidual  stoss  changes  is  monotomc  (Fig 
1.1).  I  his  signifies  that  the  matiix  mtluencevl  solelv  bv  the  Intel- 
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face  effect  is  primarily  restricted  to  a  very  small  region  in  the 
close  vicinity  of  the  whisker-matrix  interface. 

The  Volume  Mismatch  Model 

The  term  volume  mismatch  is  derived  from  the  difference  in 
volume  of  the  matrix  as  a  result  of  the  ACTE  effect  between  the 
reinforcement  and  the  matrix.  Considering  that  the  thermal  re¬ 
sidual  stress  is  in  compression  at  the  tip  of  the  whisker,  the 
material,  therefore,  may  be  regarded  as  being  composed  of  three 
distinct  zones:  the  region  of  matrix  which  is  under  tension  (de¬ 
noted  by  subscripts  ml),  the  region  of  the  matrix  which  is  under 
compression  (denoted  as  me),  and  the  whisker  zone  (denoted 
by  w').  The  combination  of  the  two  regions  in  the  matrix  ( mt  + 
me)  makes  up  the  entire  matrix  region  in  the  unit  cell. 

Although  the  local  field  in  each  zone  is  not  uniform,  since  our 
interest  is  the  average  residual  stress,  it  is  often  convenient  to 
work  with  volume  averages  of  the  non-uniformed  field  in  each 
zone  defined  by  Eq  11,  which  is  equivalent  to  Eq  3  defined  in  a 
previous  section: 

<•>*.  =  j-  fAJ  )Jy  (ID 

where  An  is  tne  area  of  each  zone,  that  is,  n  =  mt.  me.  and  w 
Considering  virtual  work  and  Gauss'  theorem,  it  can  be  shown 
that  the  following  relationship  exists: 

JJ  tT„dA  =  j  (Tlknkx,dl  (12) 

that  is, 

J  (13) 

where  nk  is  the  normal  to  the  boun.ar;  ' 

Furthermore,  it  is  common  prav-  “  to  assume  th  t  tl.c  average 
strain  ((r„>,,„)  and  avtiage  stre  ,  )  are  related  in  the  same 

way  as  in  a  loc.V  "eld  by  the  sa,  •  , ’ferial  constants  [20 1. 

For  matrix  material. 


From  the  above  simplifications,  it  is  easily  revealed  that  the 
concept  of  volume  mismatch  has  reduced  to  the  mismatch  be¬ 
tween  the  whisker  and  the  matrix  at  the  whisker  end.  Because 
all  the  quantities  we  are  interested  in  are  along  the  A\  axes,  we 
may  substitute  a  simpler  notation  for  the  element  to  be  inves¬ 
tigated,  that  is, 

(G.iXs,  =  <£/.„>;  <<To>a„  ~  (16) 

It  is  clear  that  Eqs  12  and  13  define  the  relationship  between 
average  stress  over  the  area  of  one  particular  zone  and  the  line 
average  of  the  stress  on  the  boundary  |A„|.  Theretore.  global 
equilibrium  condition  (equilibrium  between  zones)  can  be  ex¬ 
pressed  in  terms  of  volume  average  stress.  The  equilibrium  within 
each  zone  is  automatically  satisfied  •  hen  adapting  Eq  13 
Furthermore,  during  the  entire  analysis,  the  continuity  condition 
at  the  interfaces  between  zones  is  .-.atiTied  on  an  'erage  basis. 

From  the  requirement  of  the  condition  of  equilibrium  of  stresses, 
it  may  be  shown  that  the  average  residual  stress  in  the  whisker 
is  related  to  the  average  residual  stress  in  the  matrix  by  the 
following  equation: 

I  ■  VK 

<oL)  - r. —  («■;„>  ( 17) 

v  „ 

where  again.  V'„  refers  to  the  volume  fraction  of  the  whiskers 
and  1  -  V„  =  Vm,  +  Vm,  -  Vm  denotes  the  total  volume  fraction 
of  the  matrix.  The  overall  average  residual  stress  in  the  matrix 
is,  in  turn,  related  to  the  residual  stress  in  the  tensile  part  of  the 
matrix,  (it'..  '• 

(<r'J  =  F,<<0  ( IK) 

where  /•„  is  defined  by: 


(n„>A,  -  (14) 

>A,  ~  C.J  +  r  it;  '|<rrA/).,„  -  <r';,|K<r,,>.,„  -  (<rn)J 

where  C"k,  and  t.  „7)  are  the  stiffness  and  the  tangential  stiffness 
of  the  matrix,  and  f[]  and  tr1,)  are  the  average  strains  and  stresses 
on  the  yield  surfaces  in  the  strain  and  the  stress  space,  respec¬ 
tively.  and  Tresca  yield  criteria  is  used. 

Yano  et  al.  [37]  stated  that  the  lateral  stress  generated  by  the 
difference  in  Poisson's  ratio  is  small  with  a  AI  AI,  Ni  fiber-rein¬ 
forced  eutectic  composite  system  which  has  a  Poisson's  ratio  of 
v,  0.16  and  i\„  =  0.33.  Since  a  SiC'  Al  system  has  a  Poisson's 
ratio  of  vSl(  -  0.19  versus  i>AI  -  0.33.  we  shall  adopt  the  state¬ 
ment.  Therefore,  if  only  the  longitudinal  normal  stress  is  con¬ 
sidered.  Eq  14  will  reduce  to  the  following: 


(«r,,),. 


.  .  (<»„>.,„  “  <Tt„  <r,„ 

. . — nz —  z 


A„  ml.  me 


<o,,)1b  '  1T,„, 

■  <r,,„  (13) 


It  may  also  be  shown  that  the  residual  stress  in  the  compressive 
part  of  the  matrix  is  given  by: 


-  -j - ^7-  <<K„)  (20) 

Considering  the  cooling  process,  it  is  known  that  the  mismatch 
of  thermal  expansion  coefficients  is  large  enough  to  induce  plastic 
deformation  around  the  reinforcement  [.7A'| .  By  considering  equi¬ 
librium  and  compatibility,  in  an  average  sense,  it  can  be  seen 
that  the  following  equations  can  be  obtained  when  [A 7',  5  iA/  i 
<  lATj  (the  details  of  the  derivation  are  shown  in  Appendix  B 
and  in  Ref  3b)  where: 


I A  7,| 


I  AT, 


(1  4  h 


2i'„ „  /; 


i  -  i'„„  i  v  „  /  . 
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here,  |A7’,|  and  |A7',.|  represent  the  magnitude  of  the  temperature 
change  when  matrix  tensile  zone  (Aml)  and  compressive  zone 
(Amc)  become  fully  plastically  deformed,  respectively,  by  the 
thermal  stresses.  Therefore,  the  thermal  residual  stress  along  the 
longitudinal  direction  is: 


<<o 


V  F 
i  -  v; 


a% 
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(1  -  F.)  +  F. |Aa|  |A7~[ 


(22) 


When  the  applied  load  is  compressive: 


(<0  =  {a’J 


1  -  VC 


Et 


(1  -  E„)E„  +  F„E„ 

(T 


E,I= _ (iLtV 

FJE„  +  FnEm  \Ej 


(23  d) 


Equation  23  clearly  indicates  that  (<C)  is  no  longer  a  monotonic 
function.  These  equations,  therefore,  predict  a  non-monotonic 
change  in  the  residual  stresses  in  the  composite  depending  on 
whether  the  system  is  subjected  to  a  compressive  or  tensile  uni¬ 
axial  loading. 

By  considering  applied  composite  stress  (<j  J  and  the  corre¬ 
sponding  composite  strain  (f  „„  = 


This  means  that  within  a  certain  range  of  temperature  variation 
(IA7VI  £  I A  7"|  <  A 7,1 )  plastic  detormation  will  be  introduced  in 
the  matrix.  However,  until  jA 7~|  2  |A7',|,  the  matrix  will  remain 
partially  plastically  deformed.  For  ]A7'|  =  500°  C.  it  is  both 
predicted  by  Eq  2  and  experimentally  observed  by  Vogelsang  et 
al.  [39]  in  terms  of  generation  of  dislocations  that  the  matrix  is 
partially  plastically  deformed.  The  predicted  average  residual 
stress  in  the  matrix  is  given  by  Eq  22. 

It  is  of  particular  interest  to  evaluate  the  new  residual  stresses 
in  the  composite  when  a  uniaxial  loading  (it,.,)  corresponding  to 
the  total  applied  strain  (i„„)  is  applied  along  the  axis  of  the 
whisker  and  is  subsequently  removed.  Clearly,  the  applied  stress 
changes  the  residual  stress  when  plastic  deformation  is  involved. 
If  the  conditions,  as  stated  before,  are  such  that  only  portions 
of  the  matrix  have  initially  yielded  because  of  the  thermal  resid¬ 
ual  stresses,  it  may  be  shown  that  when  the  ap.  lied  load  is  tensile, 
the  overall  average  residual  stress  in  the  matrix  after  unloading 
is  given  by  (for  derivation  details,  refer  to  Appendix  B) 


<r,„  =  V  E„f,u,  (24) 

where  E,,  is  the  Young's  modulus  or  composite  tangential  Young's 
modulus  of  the  composite  depending  on  the  nature  of  the  ‘-train 
(f,„).  'hat  is.  elastic  or  plastic.  The  result  is  shown  in  Fig.  IX.  It 
is  clear,  bv  comparing  with  Fig.  4.  that  the  most  prominent 
features  predicted  by  the  FEM  are  reproduced  by  the  volume 
mismatch  model. 

Using  Eq  23.  the  average  values  ol  the  residual  stress  in  the 
matrix  were  plotted  in  Fig.  19  after  the  application  of  tensile  and 
compressive  loadings.  For  the  purpose  of  comparison,  the  results 
of  the  finite-element  analysis  described  before  have  also  been 
superimposed  on  the  same  diagram. 

-1 20 - 


and  E,,  represents  the  all-elastic  modulus  of  the  composite  in  I  K;  is  s wrwm  ,  p,,;h,  ud  l<\  ih,  i  ,,,’amr  nn\m,iu  h  m. ></.■/ 

the  absence  of  anv  plastic  thermal  residual  stresses.  «  hu  h  arc  m  ;rmw/ <wn; ni,-ni  »  nh  sh <  /  /  M  pr,  ,h,  :>•«  \hnn  >:  in  A/e  j 
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FIG  1 4  —  (  ompartson  of  average  longitudinal  residual  stress  preda  ted 
by  HIM  and  the  volume  mismatch  model.  Sole  the  similarities  of  ■ lope 
change  between  the  two  models. 


Discussion 

The  results  of  the  FEM  modelling  that  have  been  duplicated 
by  the  phenomenological  modelling  indicate  that  residual  stress 
produced  by  the  ACTE  and  cooling  is  affected  by  the  external 
load  even  at  the  initial  loading  stage.  If  the  deformation  is  a 
result  of  compressive  loading,  the  residual  stress  at  first  decreases 
and  then  increases  with  increasing  deformation,  whereas,  if  the 
deformation  is  the  result  of  tensile  loading,  the  residual  stress 
decreases  initially  at  a  rapid  rate  and  the  rate  slows  upon  further 
deformation.  If  the  tensile  deformation  is  continued,  the  residual 
stress  changes  sign  from  tensile  to  compressive. 

From  the  phenomenological  analyses,  it  is  clear  that  the  vol¬ 
ume  mismatch  model  gives  qualitative  agreement  with  the  finite- 
element  analysis  on  the  general  trend  of  the  residual  stress  change. 
Not  only  the  immediate  drop  of  the  average  residual  stress  upon 
loading  is  apparent,  but  the  slope  change  can  also  be  predicted. 
Note  that  in  the  case  of  compressive  loading,  there  is  a  sign 
change  in  the  slope  after  the  composite  has  started  yielding, 
which  was  not  predicted  by  the  interface  model  Furthermore, 
bv  comparing  Figs  4  and  IS,  the  constitutive  behavior  of  the 
composite  predicted  by  FEM  and  the  volume  mismatch  model 
are  generally  in  agreement. 

It  should  be  pointed  out  that  the  predicted  difference  in  the 
apparent  Young's  modulus  under  a  loading  scheme  of  different 
sense  is  a  direct  result  of  the  asymmetric  characteristic  of  the 
material  under  plastic  deformation.  When  an  additional  defor¬ 
mation  is  imposed  on  a  material  that  is  already  under  plastic 
deformation,  the  relationship  of  the  deformation  to  the  applied 
load  dep  i  on  whether  the  applied  load  represent  a  process 
of  loading  >  ;  unloading.  In  the  composites,  the  matrix  is  plas¬ 
tically  deformed  as  a  result  of  thermal  residual  stress  As  can  be 
seen  in  this  analysis,  the  nature  of  the  longitudinal  residual  stress 
(tensile  or  compressive)  is  distinctly  different  depending  on  the 
location  in  the  matrix  When  an  external  load  is  applied,  whether 
it  is  a  process  of  loading  or  unloading  at  different  locations  in 
the  matrix  depends  on  the  direction  of  the  external  load.  There¬ 
fore,  the  apparent  Young's  modulus  of  the  composites  is  different 
when  the  loading  direction  is  different. 


To  understand  the  underlying  mechanism  for  the  change  in 
residual  stress,  it  is  helpful  to  retrieve  the  argument  of  weighted 
residual  stress  in  the  section  of  FEM  analysis.  Since  the  ratio  of 
the  simple  average  stress  over  weighted  average  stress  is  greater 
than  I  (1.55),  this  means  that  the  magnitude  of  the  weighting 
factor  (stress  rate  bajdq)  is  greater  in  the  compressive  region. 
Physically,  it  means  that,  to  maintain  compatibility  between  the 
matrix  material  at  the  whisker  tip  and  the  rest  of  the  matrix,  the 
compressive  zone  at  the  whisker  tip  responds  faster  than  the  rest 
of  the  matrix,  that  is.  higher  stress  (strain)  rate.  This  will  result 
in  localized  plastic  deformation  before  global  yield.  Since  the 
gradient  of  plastic  deformation  in  the  material  is  mainly  respon¬ 
sible  for  the  deformation  generated  residual  stress,  the  imme¬ 
diate  change  of  the  residual  stress  upon  the  initial  external  de¬ 
formation  is  explained. 

It  is  also  interesting  to  note  that  the  current  Eshelbv  method 
does  not  predict  such  a  residual  stress  reduction  after  an  initial 
prestrain.  Withers  et  al.  \20\  in  their  recent  mean  field  investi¬ 
gation  implied  that  the  residual  stress  stays  as  a  constant  within 
a  preglobal-yield  region  and  that  there  is  an  overall  monotonic 
change  in  residual  stress  as  a  result  of  prestraining.  The  discrep¬ 
ancy  lies  in  that,  in  the  mean  field  theory,  the  yielding  criteria 
is  satisfied  in  an  average  sense.  However,  well  before  the  average 
global  yielding  criteria  is  satisfied,  localized  microyielding  has 
already  occurred,  for  example,  at  the  whisker  end.  as  pointed 
out  previously,  and  incompatibilities  between  different  parts  of 
the  matrix  which  are  responsible  for  the  generation  of  the  matrix 
residual  stress  have  already  developed.  Therefore,  residual  stress 
alteration  is  induced  before  the  global  yielding.  However,  by 
considering  only  a  global  average  field,  such  localized  micro¬ 
plastic  deformation  can  not  be  taken  into  account.  In  the  mean 
field  theory,  assumption  was  made  that  deformation  was  uniform 
in  the  matrix.  Therefore,  deformation  generated  matrix  residual 
stresses  are  due  only  to  the  incompatibility  between  matrix  and 
the  reinforcement.  The  residual  stresses  as  a  result  of  the  plastic- 
deformation  gradient  in  the  matrix  are  not  considered. 

In  the  postyield  regime,  it  is  worthwhile  to  point  ,  ut  that  in 
the  case  of  compressive  loading,  there  is  a  sign  change  in  the 
slope  (i<<r;„>/rtft„t  after  the  composite  yielding,  which  is  predicted 
by  both  FEM  and  the  volume  mismatch  composite  model.  This 
is  an  indication  that  the  interaction  between  the  tensile  and  com¬ 
pressive  zones  becomes  totally  different.  The  original  compres¬ 
sive  load  will  produce  a  geometric  compromise  between  the  ten¬ 
sile  and  compressive  zones,  and  therefore,  there  will  be  a  reduction 
in  residual  stress.  After  the  matrix  material  is  completely  de¬ 
formed  plastically  by  the  applied  compressive  load,  the  two  zones 
become  more  and  more  incompatible.  Therefore,  the  residual 
stress  increases  again.  It  is  also  important  to  point  out  that  while 
the  trend  of  the  residual  stress  change  is  in  agreement  with  both 
FEM  and  volume  mismatch  models,  the  actual  magnitude  of  the 
residual  stresses  predicted  by  the  analytical  model  is  much  smaller 
than  that  predicted  by  the  FEM.  Furthermore,  the  relative  po¬ 
sitions  between  tensile  and  compressive  curves  predicted  by  FEM 
are  not  in  agreement  with  the  volume  mismatch  model.  It  is 
speculated  that  the  absence  of  interfacia)  bonding  in  the  model 
is  the  key  reason  for  the  discrepancy  .  For  example,  the  discrep¬ 
ancy  of  relative  positions  of  the  curves  of  residual  stress  after 
tension  and  compression  between  the  FEM  and  volume  mis¬ 
match  can  be  explained  by  the  fact  that  the  effect  of  interface 
will  enhance  the  initial  residual  stress  drop  when  under  external 
tensile  load  and  oppose  such  a  drop  when  under  compression 
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Such  an  effect  can  be  verified  easily  through  Eq  10,  which  is 
directly  derived  from  shear  load  transfer.  However,  the  pattern 
of  the  residual  stress  changes  is  solely  determined  by  the  volume 
mismatch  mechanism. 

During  the  2-D  FEM  analysis,  an  assessment  was  made  that 
a  2-D  plane-strain  FEM  analysis  is  a  more  economical  way  of 
modelling  the  generation  of  the  residual  stress.  The  incompati¬ 
bility  between  matrix  material  at  the  tip  of  the  whisker  and  the 
rest  of  the  matrix  during  plastic  deformation  is  responsible  for 
the  change  in  the  residual  stress.  In  a  2-D  plane-strain  approx¬ 
imation,  this  factor  still  remains  because  of  the  existence  of  the 
tensile  and  compressive  zones  in  the  matrix.  Whereas,  in  a  more 
realistic  3-D  model,  the  only  thing  that  makes  it  closer  to  reality 
is  that  the  surface  area  of  the  reinforcement  is  significantly  in¬ 
creased  for  the  same  volume  fraction,  which  means  that  the 
mechanism  of  load  transfer  becomes  more  effective.  However, 
it  has  been  demonstrated  that  the  effect  of  load  transfer  is  only- 
restricted  to  a  small  volume  near  the  matrix-reinforcement  in¬ 
terface  in  a  whisker-reinforced  composite,  and  the  effectiveness 
of  the  load  transfer  does  not  change  the  fact  that  there  is  a  higher 
stress  concentration  at  the  tip  of  the  whisker.  Therefore,  in  a 
2-D  plane-strain  condition,  while  all  the  unnecessary  complica¬ 
tions  are  eliminated,  the  main  mechanism  of  residual  stress  gen¬ 
eration  during  plastic  deformation  is  preserved.  Considering  the 
overall  effect  of  interfacial  load  transfer  and  the  effect  of  volume 
mismatch,  it  is  suggested  that  the  overall  increase  in  the  mag¬ 
nitude  of  the  slope  in  the  change  of  residual  stress,  shown  in  Fig. 
15.  is  due  to  the  enhanced  effect  of  interfacial  load  transfer.  This 
enhanced  load  transfer  is  believed  to  he  a  result  of  the  increase 
of  matrix-whisker  interfacial  area. 


Conclusions 

In  summary: 

1 .  The  predicted  change  of  the  average  residual  stress  is  not 
monotonic  with  respect  to  the  applied  prestress  (strain),  that  is. 
regardless  of  the  sign  of  the  loading.  The  average  residual  stress 
drops  upon  the  initial  elastic  deformation,  and  the  change  of  the 
residual  stress  as  a  result  of  the  external  load  depends  on  whether 
the  loading  history  is  compressive  or  tensile. 

2.  The  results  predicted  by  FEM  arc  generally  comparable  to 
those  predicted  by  the  theoretical  model  which  is  based  on  the 
volume  mismatch  between  the  reinforcement  and  the  matrix 
This  indicates  that  the  volume  incompatibility  is  the  controlling 
factor  regarding  the  change  in  residual  stress. 

3.  The  FEM  and  theoretical  results  indicate  that  the  asym¬ 
metrical  constitutive  behavior  of  whisker-reinforced  SiC'Al  com¬ 
posites  is  primarily  due  to  the  presence  of  the  thermal  residual 
stress. 

4.  A  2-D  plane-strain  analysis  is  an  adequate  and  yet  more 
economical  approach  than  the  actual  time  consuming  3-D  FEM 
analysis  in  understanding  the  generation  of  the  residual  stress 
caused  by  external  deformation. 
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In  this  appendix,  a  modified  shear  lag  model  is  derived.  A 
more  complete  form  can  be  found  in  Ref  35. 

Considering  equilibrium,  the  follow  ing  equation  can  be  obtained: 


Assuming  that  the  aluminum  matrix  follows  the  bilinear  stress- 
strain  relationship  and  considering  displacement-strain  transfor¬ 
mation  equation,  the  displacement  can  be  related  to  interfacial 
shear  stress  in  the  following  manner: 

/»<|Auj  \AT\z  +  An  -  Au.)  -  ~  -  t.  )  ( A2) 

Combining  Eqs  A!  and  A2: 


considering  it,,,,  >  Aw,,  and  solving  Eq  A3  in  a  similar  fashion 
as  described  in  Ref  35.  it;,  can  be  obtained  as: 

/  cosh  Be  \ 

<  -  /:,,(< t:„>  -  A..:  A/  )  I - ( A4 ) 

'  cosh  13/  ' 

where  <»•;„»  is  the  averaging  matrix  thermal  residual  strain,  and 
(3  can  be  defined  as: 


Averaging  <f„  by  means  of  the  continuum  form  of  Eq  3: 


„ 


. (?) 


(?) 


(<»:.,>  Au  A/  l  (As) 


where 

■T.  4’  •  or;, ,.'(  1  l'„ )  (I  (Ah) 

Rearranging  Eqs  A5  and  Ah  and  considering  the  following: 


u;„> 


ir„  or.„>  IT 


/. 
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the  following  equations  can  be  obtained: 


APPENDIX  B 


«>  =  E„  I  1  - 


,anH(f) 


\\  E. 
1  -  V.  E„ 


«)  i'.  -  (<o(i  -  vj  -  o 

It  is  likely  that  the  matrix  is  fully  plastically  deformed  as  a 
result  of  thermal  residual  stress.  Therefore,  h  and  3  can  be  de¬ 
termined  as: 


Volume  Mismatch  Model 

In  this  Appendix,  the  volume  mismatch  model  is  derived.  A 
more  complete  form  can  be  found  in  Ref  35. 

On  satisfying  global  compatibility,  the  thermal  residual  stress 
can  be  obtained.  If  the  temperature  change  is  small  so  that  the 
matrix  remains  elastic,  the  global  compatibility  and  equilibrium 
condition  can  be  obtained  as: 


(1  _  f  j  _  F  =  |Aaj  (ATIT. 
£m  £  '  En 


(O'nu)  =  - 


The  corresponding  matrix  thermal  residual  stress  is: 


2V2 


(<T’„)  = 


V,„,F„  _ lAa||A71F. _ 

1  -  v*  1  1  -  f--  ,  F- 

£„  +  £„,  1  -  V„,  £„  1  -  V„ 


When  a  longitudinal  applied  stress  (<r  ,)  is  imposed,  the  matrix 
is  plastically  deformed  when  the  matrix  stress  is  larger  than  the 
yield  strength  of  the  matrix.  Considering  the  corresponding  ap¬ 
plied  plastic  strain  the  residual  stress  after  unloading  can 

be  determined  as: 

=  <<r;„>  +  f  ,£,„.,  -  (f, (AH) 


(a'“)  =  (o;,)  +  £„f,„;,  I  I 


£ .  'f. I  1 


— — — —  <<r'"> 

(1  -  V'J 


where  (f is  the  average  strain  resulting  from  unloading. 

Eliminating  (t,ml . ,>  from  Eqs  AH  through  All),  the  final  form 

of  the  residual  stress  is  obtained: 


<«C>  “  l',rL>  f . , I  I 


lanh(|') 

i?r 


2\  2  I  a..,,.  /,  „ 


When  the  temperature  change  is  small,  there  is  no  thermally 
induced  plasticity  in  the  matrix  When  the  temperature  change 
is  smaller  than  the  critical  temperature,  the  matrix  remains  elas¬ 
tic.  that  is. 

1 A  71  -=  A/ 

where  the  critical  temperature  change  (  AT.  )  can  be  obtained 


(ii-  i'„„  1  -  £.,  .  £ 


\  E„  T.,„ 


£,.,  £.,  1 


When  the  temperature  change  is  large  enough,  the  matrix  at 
the  tip  of  the  whisker  starts  to  yield.  The  variation  ot  the  tem¬ 
perature  change  for  such  a  condition  if  given  by: 

A/.  -  A/  •  AT 


U/;;  - 


::  -  £„)  i  -  2T„„  i  /■„  f„„  £.. 

£„,  “  1  -  V„„  £  *  1  -  V 


Considering  a  compatibility  condition  similar  to  Eq  Bl.  the 
thermal  residual  stress  can  be  obtained: 


<«r'  )  - 


<r, -  -M(l  £.,)  -  £..  Aa  A  / 

T  /•  '£»./■ 

'►If1  w _ _ _ 

’■  r-  —  .  _L=U!_L  .  L) 

£..,  I  -  l-  A  £ .  £,.  1 


I  -  l, 
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For  an  elastic  matrix,  when  an  external  load  is  applied,  the  fol¬ 
lowing  equations  can  be  constructed: 

7?  =  IT  0  -  F.)  +  7^  K 
Ect  fc* 

=  <<0(1  -  Vml)  +  <o.„)  (B6) 

Vmf.  <Q  = 

\  -  Vm  Em  £,, 

where  <crH)  and  <a„)  are  average  stress  in  the  whisker  and  matrix, 
respectively.  Solving  Eq  B6,  the  Young's  modulus  can  be  ob¬ 
tained  as: 

(1-1/  IE  £ 

=  — U - —  +  y  £  (B7) 

"  £.0  -  EJ  +  F„Em 

With  a  presence  of  thermally  induced  matrix  plasticity,  when 
an  applied  load  is  imposed  on  the  composites,  the  stress  state  in 
the  composites  can  be  expressed  as: 


For  a  tensile  load,  the  following  equations  can  be  obtained: 


1  -  Vw 


<<0  +  E„feup  Frup>(l 


<<o  = 


(1  -  F„)E,  + F„Em  "" 

~2(<>  +  (i  -  f,)E,  +  £„£„,/“"’ 


1  -  1/J2|(1  -  £„)£»  +  E.EJ  1  -  Vm,\(<rj  o.„ 

E“"’1  Ek  {  £.  V„„  J  E„  £,„ 

When  the  composites  are  under  compression,  the  internal  stress 
is  given  by: 


*  ,r'~) ' 


(<rj  =  <<0  +  OO 

(<rm>  =  <(T,'n>  +  Or",) 


where  (<0  and  Or",)  are  the  average  stresses  in  the  whisker  and 
the  matrix  as  a  result  of  applied  load,  respectively,  and  their 
value  is  different  under  different  loading  directions. 

The  total  applied  strain  can  he  related  to  the  average  whisker 
strain  and  the  strain  in  the  matrix  at  the  tip  of  the  whisker: 

£„„  -  FJrJ  +  (1  -  £„  Xi.,,.)  (B9) 

Equation  Bd  has  provided  a  relationship  between  average  strains 
from  different  locations  in  the  matrix. 

To  consider  the  stress  gradient,  the  matrix  average  stresses  in 
the  initial  compressive  zone  ((it,,,.))  and  tensile  zone  should 

be  considered  separately.  Using  the  equilibrium  condition,  the 
following  equations  can  he  obtained: 


V,n,  V',„ 

0r,„>  -  -j - —  <ir„„)  +  -j - —  <<r,„.  > 

<r.„  ~  (I  -  V’„|(<T,„)  + 


Because  of  the  thermally  induced  matrix  plasticity,  the  re¬ 
sponse  of  the  composites  to  a  uniaxial  tensile  load  is  different 
from  a  compressive  load  Considering  Eqs  Bd  and  Bln.  and 
following  the  development  of  the  matrix  plasticity,  the  response 
of  the  composites  to  the  external  applied  strain  can  be  derived 


<<C,>  ~  (1  -  FJ£„  +  £„£„,/"" 

After  the  external  applied  stress  is  removed  from  the  com¬ 
posites.  deformation-induced  residua!  stresses  in  the  matrix  can 
then  he  determined  Considering  Eq  B8.  the  residual  stress  state 
can  be  obtained: 

(<<C>  -  (<r'„>  *  Ut:>  +  BM 

|  Or;;;)  -  Or;;,)  +■  ut;""""1/ 

where,  by  considering  liqs  15  and  BK  through  Bid.  and  paying 
special  attention  to  the  stress  state  at  each  zone,  that  is.  Or„„>. 
0r„>.  and  0t„„>.  the  matrix  stress  as  a  result  of  unloading 
can  be  determined  as: 


l  -  U„  (1  -  £.,)£„  *  £„£„,,.  £  , 

Substituting  Eqs  B5.  B1 1  (or  BI2).  and  B14  into  Eq  B 1 3 .  and 
with  a  consideration  of  Eqs  Bd  and  Bid.  the  matrix  residual 
stress  can  be  generated 
For  compressive  load: 


dr’'')  -  Or' ) 


cj  (1  -  £,)£..  -  £..£.,„ 

£„£,„  <r„l 

(1  £..)£..  £„/  ...  £.,  ' 


For  tensile  load: 


„ _ £„£„, 
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Fracture  toughness  of  discontinuous  metal  matrix  composites 

R.J.  Arsenault 

Metallurgical  Materials  Laboratory,  Materials  and  Nuclear  Engineering  Department, 
University  of  Maryland,  College  Park,  MD  20742-2115 

Abstract 

The  various  theories  that  have  been  proposed  to  account  for  fracture  toughness 
of  discontinuous  metal  matrix  composites,  or  the  lack  there  of,  are  examined.  It  has 
been  observed  that  the  fracture  toughness  is  independent  of  particle  size  and  the 
fracture  proceeds  through  the  matrix.  A  significant  factor  which  has  not  been 
considered  in  the  past  is  that  the  matrix  is  in  a  severely  cold  worked  condition.  The  cold 
worked  state  is  the  major  factor  contributing  to  the  low  fracture  toughness  of 
discontinuous  metal  matrix  composites. 

1.  Introduction 

The  fracture  toughness  of  discontinuous  metal  matrix  composites  (DMMC)  is 
always  less  than  the  fracture  toughness  of  the  matrix,  in  the  same  heat  treatment  or 
annealed  condition.  If  a  specific  case  is  considered,  a  20  V%  SiCp/6061  composite  aged 
to  the  T-6  condition  has  a  fracture  toughness  (Kjc)  of  -  10  MPa-  mfi,  where  as  0  V% 
6061  may  have  a  Kj£  of  >  40  MPa  -  m*.  There  may  be  obvious  reasons  why  the  Kjc  of 
the  composite  is  less  than  that  of  matrix  alloy,  but  there  are  still  numerous  questions 
associated  with  the  fracture  characteristics  of  DMMC. 

During  the  past  several  decades  there  have  been  a  considerable  number  of 
investigations  of  the  fracture  characteristics  of  two  phase  alloys  where  the  second  phase 
is  a  discrete  precipitate  or  particulate[l-18).  The  precipitate  particle  is  usually  assumed 
not  to  plastically  deform  and  the  matrix  is  assumed  to  be  ductile.  DMMC  meet  the 
criterion  of  these  assumptions  very  well.  The  review  by  Schwalbe  [13]  and  more  recent 
publications  by  Bates  [15],  Gerberich  [16]  and  Firrao  and  Roberti  [17]  are  very 


informative  on  this  subject.  In  ductile  fracture,  characterized  by  void  nucleating  and 
growth  (VNG),  the  spacing  between  the  void  nucleating  particles  is  generally 
considered  to  be  a  critical  microstructural  parameter  and  taken  together  with  the 
tensile  properties  controls  the  toughness  of  a  given  material  [11-17],  This  can  be 
represented  by  the  following  simplified  expression  [17] 

KIC  =  [ao/fErfN)]^*  (1) 

where  Kjc  is  the  plane  strain  fracture  toughness,  a  is  a  numerical  coefficient,  oy  is 
the  yield  stress,  ef  is  the  maximum  strain  acting  at  the  crack  tip,  E  is  the  Young’s 
modulus,  f(N)  is  some  function  of  the  strain-hardening  exponent  and  s  is  the  average 
inclusion  spacing  in  the  matrix.  The  experimental  data  collected  by  various 
investigators  are  usually  given  in  terms  of  the  volume  fraction  of  the  second  phase 
particles  (V)  or  less  often  in  terms  of  s  (see  Fig.l) 

Therefore,  it  appears  that  the  fracture  toughness  of  SiC/Al  composites  can  be 
improved  when  the  size  of  the  SiC  particles  is  increased,  providing  that  at  least  one  of 
the  two  following  assumptions  is  met: 

(1) .  voids  are  nucleated  by  SiC  particles  and/or 

(2) .  the  response  of  the  matrix  to  change  of  s  remains  the  same  regardless  of 
whether  this  change  is  caused  by  varying  the  volume  fraction  or  the  size  of  the  SiC 
particles. 

Crack  initiation  and  propagation,  i.e.  the  separation  mechanism,  in  SiC/Al 
composites  has  been  investigated  to  some  extent  over  the  last  few  years.  Essentially, 
there  are  three  different  questions  that  have  been  raised  concerning  the  crack  initiation 
and  propagation.  These  are:  (1)  Does  the  fracture  of  the  SiC  particles  contribute  to  the 
separation  of  the  composite,  i.e.  are  they  an  integral  part  of  the  fracture  mode?  (2) 
Does  void  nucleation  at  the  SiC  whisker  or  particle  and  growth  contribute  to  the 
fracture  process?  (3)  What  is  the  significance  of  cracked  tip  opening  displacement? 

The  fracture  of  SiC  whiskers  or  particles  can  conceivably  initiate  fracture  of  the 
composite.  If  the  particle  or  whisker  were  large  enough,  then  a  crack  of  critical  length 
could  form  and  if  the  dynamic  Kjc  of  matrix  were  low  enough,  the  crack  would 
continue  to  propagate  into  the  matrix  leading  to  failure  of  the  composite  [19].  Embury 
et  al  [20]  have  also  developed  a  damage  initiation  model,  i.e.  cracking  of  SiC  particles, 
as  a  fracture  mechanism  of  SiC/Al  composites. 
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Fig.l  Fracture  toughness  as  a  function  of  "process"  zone  size.  When  VNG 
mechanism  is  active  the  process  zone  size  is  equated  with  the  average  spacing, 
s,between  void  nucleating  particles(Gerberich,  Ref.(16]).  Various  forms  of 
equation(l)  are  employed  by  Gerberich  to  predict  the  fracture  toughness 
values  (solid  and  dashed  lines)  and  compare  them  with  the  experimental  data 
reported  by  Schwalbe  and  Spitzig. 

If  the  fracture  of  SiC  could  initiate  fracture,  then  the  fracture  of  SiC  should  aid 
in  the  propagation  of  the  crack.  It  has  been  shown  that  large  intermetallic  inclusions 
can  act  as  a  fracture  initiation  site  in  SiC/Al  composites  [_1].  However,  in  an 
experiment  to  measure  the  bond  strength  between  SiC  and  Al,  a  tensile  sample  of 
99.99%  Al  was  fabricated  which  contained  a  3  mm  diameter  sphere  of  SiC  in  the  center 
of  the  gage  section.  The  gage  section  was  25  mm  long  and  had  a  diameter  of  6  mm. 


When  the  sample  was  tested,  shortly  after  macro-yielding,  the  SiC  sphere  fractured 
perpendicular  to  the  tensile  axis  with  a  small  drop  in  load.  The  test  was  continued  and 
the  sample  finally  failed,  but  the  sample  did  NOT  fail  by  extension  of  the  crack  in  the 
SiC  into  the  A1  matrix,  even  though  there  was  a  25%  reduction  in  the  cross  sectional 
area  as  a  result  of  the  fractured  SiC  sphere.  The  sample  failed  at  another  place  in  the 
gage  length  [22]. 

Another  factor  which  should  be  examined  in  a  consideration  of  the  fracture 
toughness  of  DMMC  is  difference  in  the  microstructure  of  the  matrix  alloy  (  in  the 
composite  )  vs  the  microstructure  of  non-reinforced  matrix  alloy.  The  microstructure  of 
the  composite  matrix  is  almost  the  same  as  a  cold  worked  microstructure.  The 
composite  matrix  microstructure  has  a  higher  dislocation  density  than  the  cold  worked 
matrix  alloy.  In  the  specific  case  of  20  V%  SiCp/1 100  composite  (5pim  SiC  particles)  the 
dislocation  density  is  higher  than  that  of  1100  A1  cold  rolled  (rolled  at  room 
temperature)  95  %  reduction  in  thickness.  Therefore,  it  is  probably  only  meaningful  to 
compare  the  fracmm  toughness  of  DMMC  with  cold  worked  matrix  alloy  [23]. 

In  the  remainder  of  this  paper,  ,  a  discussion  will  be  undertaken  of  the  three 
major  points  raised  in  the  introduction. 


2.  Particle  Size  Dependance 

In  the  investigation  by  Flom  and  Arsenault  [24],  compact  tension  samples  (CTS) 
were  tested  using  the  single  specimen  J-integral  test  method  per  ASTM  standard  E  813. 
Also,  the  energy  separation  technique  (EST)  was  utilized  as  an  additional  tool  in  the 
load-unload  records  analysis  [25].  The  EST  enables  one  to  obtain  values  of  plastic,  I, 
and  elastic,  G,  contributions  to  the  total  value  of  J-integral.  It  implies  that  the  area 
under  the  load  vs  displacement  curve  which  corresponds  to  the  work  done  by  external 
force  can  be  separated  into  the  stored  elastic  strain  (potential)  energy,  Us,  the  elastic 
energy,  Ug,  released  during  crack  extension,  and  plastic  energy,  Up,  dissipated  during 
the  crack  extension,  as  schematically  shown  in  Fig.2.  A  more  detailed  description  of 
the  EST  can  be  found  elsewhere  [25].  The  rate  of  plastic  energy  dissipation,  I  = 
l/BndUp/da,  and  the  elastic  energy  release  rate,  G  =  l/BndU£/da,  represent  the 
plastic  and  elastic  parts  of  the  J-integral,  i.e.  J  =  I  +  G,  where  Bn  is  toe  thickness  of 
the  CTS  between  the  side  groove  and  a  is  the  crack  length.  Crack  initiation  fracture 
toughness  can  be  determined  as  [25] 


GEcy. 


KIC=  <77^> 


(2) 


where  Ec  is  the  composite  Young’s  modulus  and  v  is  the  Poisson  ratio  ( v  =  0.31).  Also 
a  direct  evaluation  of  Kjq  is  possible  by  taking  maximum  load  Pm  from  the 
load-unload  record  and  substituting  into  the  expression[26] 


Kiq 


P  mf(a/w) 
Bn(w)^ 


O) 


where  W  is  the  width  of  the  CTS  as  shown  in  Fig.5.4(a)  and  values  of  f(a/w)  are  readily 
available.  Crack  growth  fracture  toughness  is  evaluated  by  a  dimensionless  tearing 
modulus,  T,  which  is  equal  to  [27] 


(4) 


where  oy  is  the  composite  yield  stress  and  dJ/da  is  the  slope  of  the  stable  crack 
extension  portion  of  J  vs  a  plot  constructed  in  accordance  with  ASTM  E813.  Crack 
extension  was  determined  by  using  the  unloading  compliance  technique  [27],  In  order 
to  verify  the  calculated  crack  extension  values,  the  tested  CTS  were  exposed  to  elevated 
temperature  and  then  fractured  (heat  tinting  method).  Experimentally  measured  crack 
lengths  were  found  to  be  within  5-10%  from  the  values  calculated  by  a  compliance 
technique. 
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Fig.2  Schematic  of  the  energy  separation  principle  between  two  successive 
unloadings. 


Crack  initiation  fracture  toughness  measured  as  Kjq  and  Kjq  is  plotted  as  a 
function  of  the  average  SiC  particle  size  and  shown  in  Fig.3.  The  fact  that  both  Kjq 
and  Kjq  show  the  same  trend,  i.e.  no  dependence  on  the  size  of  the  SiC  particles, 
increases  the  confidence  in  the  results  obtained  and  supports  the  energy  separation 
method  as  a  new  and  powerful  tool. 


Fig.3  Crack  initiation  fracture  toughness  of  SiC/Al  composites  measured  as  Kjq  as  a 
function  of  average  SiC  particle  size. 


The  numerical  difference  (Kjc  *  18  MPa  -  m/i  and  Kjq  *  23  MPa  -  m^)  between 
the  values  of  Kjc  and  Kjq  can  be  explained  as  follows.  A  finite  notch  with  a  root 
radius  of  *  150jim  was  machined  in  CTS  in  order  to  start  the  crack  (in  contrast  with 
ASTM  E  399  requirement  of  a  fatigue  precracking  to  form  a  geometrically  sharp 
crack).  Thus  it  is  reasonable  to  related  to  the  microstructure.  Substituting  values  for 
Kic(r)  =  23  MPa-  m/‘,  Kjc  =  18  MPa  -  m^  and  r  ~  150  jim,  we  obtain  c  =  20  um. 
This  value  is  within  the  range  of  the  microcracking  ahead  of  the  crack  tip  in  a  2.4  /jm 
average  SiC  particle  size  composite. 

Initiation  fracture  toughness  of  a  250  /im  SiC/Al  composite  is  almost  a  factor  of  2 
less  than  the  rest  of  the  tested  composites.  This  is  apparently  due  to  the  premature 
cracking  of  250  SiC  particles. 

Crack  growth  fracture  toughness  defined  as  tearing  modulus,  T  versus  the  average 
size  of  SiC  particles  is  plotted  in  Fig.4.  The  increase  of  the  crack  growth  toughness  with 
an  increase  of  the  SiC  particle  size  means  that  more  energy  is  dissipated  during  crack 
extension  in  the  composite  with  a  larger  size  of  SiC  particles. 


Fig.4  Crack  growth  fracture  toughness  of  SiC/Al  composites  measured  as  T  as  a 
function  of  average  particle  size  at  a  constant  volume  fraction  of  20%. 


The  above  experimental  results  indicate  that  crack  initiation  fracture  toughness 
Kjc  is  independent  of  the  SiC  particle  size  (within  the  range  of  SiC  sizes  tested).  In 
order  to  understand  this  result  it  is  instructive  to  consider  local  energy  dissipation 
mechanisms  associated  with  the  crack  extension  in  SiC/Al  composite.  We  are  using  the 
term  "local"  to  emphasize  the  area  in  the  vicinity  of  the  crack  tip.  The  plastic 
deformation  in  the  bulk  of  the  SiC/Al  specimen  is  not  included  since  it  is  separated  out 
by  using  the  energy  separation  technique  (EST). 

The  elastic  energy  release  rate  G  consists  of  the  several  terms  that  account  for 
energy  dissipation  .  These  terms  can  be  identified  as 

G  =  7s  +  7v  +  7  PLzone  +  7CTS  (6) 

where  s  is  the  energy  consumed  in  creation  of  the  new  surfaces,  7V  is  the  energy 
consumed  in  the  formation  of  the  void  sheet,  7 PLzone  *s  l^e  energy  dissipated  into 
plastic  deformation  around  the  crack  tip  and  7crs  rePresents  crack  tip  shielding 
effects  described  by  Ritchie  and  Yu  [28]. 

It  can  be  shown,  that  most  of  the  energy  is  absorbed  within  the  plastic  zone 
around  the  crack  tip  (7  PLzone)-  order  to  estimate  this  energy  one  has  to  know  the 
distribution  of  the  plastic  strain  in  a  direction  perpendicular  to  the  direction  of  crack 
propagation  [29].  The  data  indicates  that  in  the  region  immediately  adjacent  to  the 
crack  the  plastic  strain  is  large,  as  evidenced  by  very  large  dislocation  density  [30].  Also 
it  was  observed  that  the  dislocation  density  is  the  same  in  the  immediate  region  of 
fracture  for  two  different  volume  fractions  of  SiC  reinforcements  [30]  (SiC  particle  size 
was  the  same).  Based  on  these  results  we  can  assume  that  as  the  particle  size  increases, 
i.e.  an  interparticle  spacing  increases  (V  =  const.),  the  density  of  dislocation  in  the 
immediate  region  of  the  crack  remains  the  same.  Therefore  7  PLzone  *s  independent  of 
particle  size. 

The  increase  in  crack  growth  fracture  toughness  (measured  as  I  and  T)  with  an 
increase  of  the  SiC  particle  size  can  be  treated  on  the  basis  of  the  relationship  between 
the  dislocation  density  and  the  spacing  of  the  SiC  particles  demonstrated  by  Arsenault 
and  Shi  [31].  From  their  work  one  can  see  that  the  increase  in  the  tearing  modulus  with 
an  increase  in  particle  size  is  the  result  of  the  increased  plastic  deformation  that  can  be 
accommodated  in  the  matrix  as  the  interparticle  spacing  increases. 

The  data  discussed  above  can  not  support  a  mechanism  based  on  VNG  at  the 
interface  of  SiC  and  the  matrix  [18,  32]  for  these  mechanism  predict  and  increase  in 
Kjc  with  an  increase  in  particle  size  at  a  constant  volume  fraction.  Also,  the  damage 


mechanism  of  Embury  et  al.  [20]  which  predicts  on  increase  in  particle  cracking  with  an 
increase  in  particle  size.  This  would  suggest  as  the  particle  size  increases  the  fracture 
toughness  should  decrease.  This  indeed  happens  when  very  larger  particles  are 
involved  i.e.  250  pm  particle  size.  However,  in  the  particle  size  range  of  0.5-20  ^m  there 
is  no  decrease  in  Kjq  which  would  indicate  that  particle  cracking  is  playing  a  very  small 
role  in  the  fracture  process. 

3.  Separation  Mode 

As  discussed  in  the  introduction,  the  mode  of  crack  propagation  is  area  of  some 
controversy.  An  investigation  was  undertaken  to  determine  how  the  crack  propagates 
[33].  Precracked,  polished  compact  tension  samples  (CTS)  loaded  in  a  mini-tensile 
stage  within  the  scanning  electron  microscope.  The  stage  was  capable  of  incrementally 
loading  the  CTS,  which  allowed  examination  of  the  crack  tip  under  load. 

The  crack  tip  region  was  then  photographed,  produce  a  map  for  the  purpose  of 
determining  where  and  which  of  the  SiC  particles  were  fractured.  Then,  the  crack  was 
advanced  by  loading  the  sample,  and  then  rephotographed. 

Composite  samples  containing  various  size  of  SiC  particulates  were  investigated; 
however,  only  a  few  of  them  will  be  considered  here. 

The  specific  case  to  be  discussed  in  detail  is  a  composite  containing  20  V%  8  pm 
SiC  particles.  After  the  crack  had  been  propagated  a  considerable  distance,  the  crack 
has  stopped  in  front  of  a  rather  large  particle  (Fig.5).  Continued  loading  of  the  sample 
results  in  no  apparent  change  in  the  crack,  except  that  the  crack  tip  has  opened  further. 
The  crack  opening  is  extremely  small;  even  at  a  distance  of  5  to  10  pm  behind  the  crack 
tip,  the  crack  opening  displacements  is  <  1  pm  and  probably  <  0.1  pm.  Also,  the  large 
particle  at  the  tip  of  the  crack  did  not  fracture,  the  crack  proceeded  around  the  large 
particle. 

Upon  further  loading  of  the  sample,  as  shown  in  Fig.6,  the  crack  tip  has  stopped 
somewhere  in  front  of  a  particle  which  apparently  is  cracked  in  three  different 
locations.  The  crack  proceeds,  as  shown  in  Fig.7,  through  one  of  preexisting  cracks  in 
the  cracked  sample.  Figures  8  and  9  are  continuations  of  the  crack. 

From  a  consideration  of  the  above  experimental  results  it  is  obvious  that  there  are 
differences  with  some  of  the  previously  generated  data  on  the  fracture  characteristics  of 
SiC/Al  composites,  but  these  differences  can  be  readily  accounted  for.  Also,  a  clearer 
explanation  of  some  of  the  previous  data  is  now  possible. 


Fig.5  A  scanning  electron  microscope  micrograph  of  a  propagation  of  a  crack  of  a  20 
V%  8  nm  SiC  particulate/1100  A1  matrix  composite. 

The  previous  data  from  matching  fracture  surfaces  and  unloaded  propagated 
cracks  supported  the  concept  that  fracture  of  SiC  particles  in  front  of  the  advancing 
crack  played  a  significant  role  in  the  fracture  process  and  can  be  readily  explained  or 
accounted  for.  The  above  mentioned  previous  data  is  based  on  the  examination  of 
postmortem  samples.  The  basic  problem  associated  with  examining  postmortem 
samples  is  that  it  is  impossible  to  know  whether  the  SiC  particle  or  whisker  was  broken 
due  to  processing  prior  to  the  actual  fracture  test.  In  the  present  investigation  there 
was  no  case  where  it  could  be  shown  that  the  advancing  crack  resulted  in  the  fracture  of 
SiC  particles. 

It  was  clearly  observed  that  as  the  crack  advances,  it  has  a  distinct  preference  to 
proceed  through  existing  cracked  SiC  particles.  This  would  be  entirely  logical,  since 
with  the  crack  advancing  through  existing  cracked  SiC  particles  there  is  no  need  to 
create  two  new  surfaces  since  they  already  exist.  Therefore,  matching  fracture  surfaces 
would  have  a  matching  "fractured"  SiC  particle.  Also,  the  observation  by  Flom  and 
Arsenault  (24]  that  the  number  of  SiC  particles  along  the  fracture  path  is  greater  than  a 
random  line  through  the  composites  is  simply  a  consequence  that  the  crack  would 
rather  propagate  through  an  existing  cracked  (fractured)  SiC  particle. 


A  scanning  electron  microscope  micrograph  of  a  propagation  of  a  crack  of  a  20 
V%  8  ;im  SiC  particulate/1 100  matrix  composite. 


Fig.8  A  scanning  electron  microscope  micrograph  of  a  propagation  of  a  crack  of  a  20 
V%  8  ^m  SiC  particulate/1 100  A1  matrix  composite. 


Fig.9  A  scanning  electron  microscope  micrograph  of  a  propagation  of  a  crack  of  a  20 
V%  8  ^m  SiC  particulate/1 100  A1  matrix  composite. 


In  terms  of  the  sharpness  of  the  crack  tip  in  the  composites,  i.e.,  lack  of  crack  tip 
blunting,  this  can  be  traced  back  to  the  lack  of  plastic  deformation  which  does  not  occur 
in  these  composites.  In  order  to  have  crack  tip  blunting,  it  is  necessary  to  have  in-plane 
plastic  deformation.  If  plastic  deformation  in  the  crack  tip  region  (which  can  occur)  is 
very  small,  there  will  be  no  blunting.  Considering  a  very  simple  case,  if  the  crack  was 
initially  very  sharp  and  then  blunted  to  a  crack  tip  opening  displacement  of  2  ^m,  as  in 
the  case  of  the  cold  rolled  6061  A1  alloy,  the  effective  plastic  strain  required  would  be 
15  -  30%  as  obtained  by  McMeeking  [34]  in  an  FEM  investigation  of  crack  tip  blunting. 
It  may  be  argued  that  the  analysis  of  McMeeking  is  not  valid  since  his  analysis  is  for  a 
stationary  crack.  The  cracks  in  the  SiC/Al  composite  propagate  in  a  discontinuous 
manner.  The  major  component  of  this  effective  plastic  strain  is  a  tensile  strain 
perpendicular  to  the  crack  plane.  However,  Davidson  [35]  believes  that  the  tensile 
strain  is  very  small  in  SiC/Al  composites.  The  tensile  strain  is  small  due  to  the  elastic 
constraint  of  the  high  modulus  (in  comparison  to  Al)  SiC  particulate  present  in  the 
composite,  and  the  fact  that  the  matrix  within  the  composite  has  a  very  high  dislocation 
density,  i.e.,  equivalent  to  a  highly  cold  worked  condition. 

Even  though  the  6061  Al  alloy  was  cold  rolled  75%,  it  is  still  possible  to  have 
significant  plastic  strain  occurring  at  the  crack  tip  in  the  plane-stress  condition.  Prior  to 
slicing  the  sample  in  half,  the  crack  propagates  in  a  nearly  plain-strain  condition.  As  a 
result  of  the  greater  plastic  strain,  there  is  significant  crack  tip  opening  (several  orders 
of  magnitude)  in  comparison  to  the  SiC/Al  composites. 

The  lack  of  void  nucleation  and  growth  at  the  SiC  matrix  interface  is  more  than 
likely  due  to  a  very  high  bond  strength  [36]  between  SiC  and  Al. 

4.  Matrix  Microstructure  Effect 

A  general  observation  is  that  fracture  toughness  of  a  metal  or  alloy  decrease  with 
increase  in  cold  work.  In  other  words,  the  fracture  toughness  is  generally  greater  than 
the  annealed  condition  as  compared  to  the  cold  worked  condition. 

The  microstructure  of  the  matrix  of  an  annealed  or  precipitation  hardened  SiC/Al 
composite  is  similar  to  matrix  cold  worked  90  %.  This  would  suggest  that  if  the  matrix 
material  was  cold  worked  90  %,  then  the  fracture  toughness  of  the  matrix  material 
would  be  comparable  to  that  of  the  20  V%  SiCp/Al  alloy  composite. 

A  preliminary  investigation  showed  (Table  1)  that  cold  rolling  6061  Al  alloy  to 
69%  resulted  in  nearly  a  factor  of  two  reduction  in  Kj£,  and  in  the  case  of  139%  cold 
rolling  there  is  further  reduction  in  Kjq.  At  this  level  the  Kjq  value  of  the  cold  rolled 
matrix  alloy  is  comparable  to  the  20  V%  SiCp/6061  Al  alloy  composite. 


Table  1 

Fracture  toughness  vs  cold  work 


Material 

Percentage  of 
cold  work 

KlC 

MPa.m^ 

6061  A1T6 

0 

43 

6061  A1 T6 

15 

31 

6061  A1  T6 

69 

27 

6061  A1 T6 

139 

-20 

In  another  investigation  [37]  in  which  correlation  between  Kjc  and  yield  stress 
was  determined,  it  was  shown  that  as  the  yield  stress  increased  (due  to  cold  work)  the 
fracture  toughness  decreased.  At  a  yield  stress  (of  the  cold  worked  matrix  alloy) 
corresponding  to  yield  stress  of  20  V%  SiCp/6061  A1  alloy  composite,  the  cold  work 
6061  A1  alloy  has  the  same  Kjc  value  as  the  20  V%  SiCp/6061  A1  alloy  composite. 

5.  Conclusions 

From  a  consideration  of  the  data  discussed  above,  it  possible  to  arrive  at  several 
conclusions. 

1.  An  increase  of  SiC  particle  size  did  not  improve  Kjc  fracture  toughness  of 
SiC/Al  composite.  Kjc  was  independent  of  SiC  particle  size  up  to  20  pm  average  size. 

2.  The  crack  growth  toughness,  i.e.  the  tearing  modulus  and  plastic  work, 
increased  as  the  SiC  particle  size  increased. 

3.  No  evidence  of  void  nucleation  at  SiC  particles  was  found. 

4.  SiC/Al  composite  is  a  more  complicated  system  than  it  appears:  (1)  the  sites  of 
void  nucleation  were  not  well  defined,  and  (2)  the  two  alternative  ways  of  increasing 
interparticle  spacing  did  not  produce  the  same  results: 

•  If  the  particle  size  was  constant,  and  the  volume  fraction  was  increased  and 
spacing  decreased,  then  Kjc  decreased. 

•  If  the  particle  size  and  spacing  were  increased,  and  volume  fraction  was  held 
constant,  then  Kjc  remained  constant. 


5.  The  crack  appears  to  go  around  intact  SiC  particles  leaving  a  distinct  outer 
coati  on  the  SiC  particles  which  is  in  agreement  with  in  situ  TEM  investigations  of 
crack  propagation  in  this  particular  material. 

6.  The  observance  of  cracking  of  the  SiC  particles  in  front  of  the  crack  tip  due  to 
the  advance  of  the  crack  tip  is  a  rare  event.  If  there  are  preexisting  fractured  SiC 
particles,  the  crack  proceeds  through  the  precracked  SiC  particles. 

7.  The  crack  tip  opening  displacement  at  the  crack  tip  and  5-15  /im  behind  the 
crack  tip  as  it  initially  propagates  in  the  composites  containing  2.4-20  vm  size  SiC 
particles  was  extremely  small,  i.e.,  less  than  0.1 

8.  An  important  factor  where  considering  the  fracture  toughness  of  SiC/Al  composites 
is  that  matrix  has  a  high  dislocation  density  and  a  small  subgrain  size,  even  where  the 
composite  is  in  the  annealed  condition.  Therefore  the  matrix  of  the  composite  can  be 
equated  to  a  cold  worked  condition.  Cold  working  always  reduce  the  fracture 
toughness. 
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Influence  of  the  Thermally  Induced  Plasticity  and 
Residual  Stresses  on  the  Deformation  of  SiC/Al 

Composites 

N.  Shi  and  R.  J.  Arsenault 


ABSTRACT 

The  Finite  Element  Method  (FEM)  was  employed  to  investigate  the  matrix  plastic  flow 
in  a  whisker  reinforced  SiC/Al  composite  under  external  tensile  load.  It  was  found  that 
the  plastic  zone  induced  by  the  plastic  relaxation  of  thermal  stresses  expands  under  the 
external  tensile  load.  The  overall  matrix  plastic  flow  was  characterized  by  the  expansion 
and  interconnection  of  the  plastic  zones  around  whiskers.  This  process  can  be  divided 
into  several  characteristic  stages,  and  related  to  the  global  stress-strain  relationship.  It 
was  also  found  that  composite  asymmetric  constitutive  response  to  external  tensile  and 
compressive  loads  as  well  as  the  composite  asymmetric  Bauschinger  behavior  was  due  to 
the  thermally  induced  plasticity  and  residual  stresses. 


1  INTRODUCTION 

Arsenault  et  al.  proposed  that  a  high  dislocation  density  due  to  the  difference  in 
coefficients  of  thermal  expansion  ( ACTE)  should  account  for  the  strengthening  of  com¬ 
posites  [1]  Miller  et  al.  suggested  that  the  strengthening  of  SiC/Al  composites  may 
result  from  the  Orowan  mechanism  |2].  However,  it  is  still  unclear  as  how  the  thermally 
induced  plasticity  affects  the  evolution  of  the  matrix  plastic  flow  which  is  important  to 
the  understanding  of  the  composite  strengthening. 

When  SiC/Al  composites  are  subjected  to  either  a  tensile  or  compressive  load,  their 
constitutive  behavior  is  distinctly  different.  It  was  generally  found  that  the  compressive 
yield  strength  of  these  composites  was  larger  than  that  of  the  tensile  yield  strength, 
while  the  apparent  Young’s  modulus  for  compression  was  smaller  than  that  for  tension 
[3].  If  a  material  is  initially  plastically  deformed  in  tension,  then  reverse  deformed  in 
compression,  a  difference  in  stress  is  required  to  initiate  plastic  deformation  in  the  reverse 
cycle.  Usually,  the  magnitude  of  the  yield  stress  in  the  reverse  cycle  is  less  than  the  flow 
stress  in  the  forward  cycle.  This  reduction  in  stress  (or  the  strain  necessary  to  reach 
the  previous  level  of  flow  stress)  is  defined  as  the  Bauschinger  Effect.  In  the  case  of 
monolithic  Al,  the  direction  of  initial  deformation  (i.e.,  tensile  or  compressive)  does  not 
have  any  influence  on  the  Bauschinger  Effect.  The  SiC/Al  composites,  however,  exhibit 
a  remarkable  characteristic,  i.e.,  the  flow  stress  drop  (or  alternatively,  the  Bauschinger 
strain)  is  larger  when  the  composite  deforms  in  a  compression-tension  sequence  than  vise 
versa  [3]. 

In  this  investigation,  FEM  modeling  was  performed  on  the  evolution  of  matrix  plastic 
flow  of  a  SiC  whisker  (SiCw)  reinforced  annealed  6061  Al  by  monitoring  the  load-induced 
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changes  of  the  plastic  zone  generated  by  ACTE,  and  then  a  correlation  was  determined 
between  the  matrix  microplastic  flow  and  the  global  composite  tensile  stress-strain  curve. 
Based  on  the  knowledge  of  the  evolution  of  the  matrix  plasticity,  a  fundamental  quest 
leading  to  an  understanding  of  the  asymmetric  composite  constitutive  behavior  was  ini¬ 
tiated,  especially  the  asymmetric  apparent  Young’s  modulus  and  the  Bauschinger  Effect. 


2  FEM  ANALYSIS 

In  this  investigation,  a  2-D  plane  strain  FEM  thermo-elasto-plastic  analysis  was  per¬ 
formed,  in  which  the  Von  Mises  yield  criteria  and  the  incremental  flow  rule  were  employed 
[4].  The  composite  was  assumed  to  be  a  2-D  infinite  periodic  array  of  hexagonally  dis¬ 
tributed  SiC  whiskers  of  20V%  embedded  in  the  matrix  (or  the  staggered  array),  as 
shown  in  Fig.  1.  An  aspect  ratio  of  four  was  selected.  This  value  is  a  typical  average 
from  experimental  observations  [5,  6].  Due  to  symmetry  conditions,  a  unit  cell,  as  shown 
in  Fig.  1,  may  be  selected.  The  details  of  the  boundary  conditions  along  the  unit  cell 
boundary  were  published  elsewhere  [7]. 

To  account  for  the  thermally  generated  plasticity  due  to  cooling  from  the  annealing 
temperature,  a  temperature  change,  AT,  was  applied  to  the  FEM  unit  cell.  The  materials 
properties  were  obtained  from  Ref.  [8,  9]  During  modeling,  neither  thermal  gradient  nor 
matrix  creep  was  considered.  Due  to  a  large  ACTE,  thermally  induced  matrix  plastic 
flow  from  AT  =  480°  C  is  likely  to  spread  over  the  entire  matrix.  To  better  understand 
the  evolution  of  the  matrix  plasticity,  two  temperature  changes  of  AT  =  30  and  480°C 
were  employed,  respectively. 

It  has  been  reported  previously  that  the  matrix  longitudinal  thermal  residual  stress  is 
compressive  at  the  tip  of  the  whisker,  and  tensile  in  the  rest  of  the  matrix  [10].  Therefore, 
in  the  current  analysis,  the  matrix  is  categorized  into  two  different  zones:  the  initial 
compressive  and  tensile  zones  at  the  tip  of  and  along  the  longitudinal  sides  of  the  whiskers, 
respectively.  The  schematics  of  this  partition  is  shown  in  Fig.  1. 


3  RESULTS  AND  DISCUSSIONS 


3.1  THE  BEHAVIOR  OF  THE  PLASTIC  ZONE 

Following  the  temperature  change  of  AT  =  30°C,  longitudinal  loading  is  applied  to 
AB  (in  Fig.  1)  incrementally.  Fig.  2  shows  the  behavior  of  the  thermally  induced  matrix 
plastic  zone  in  a  20V%  composite  under  external  tensile  loading.  In  this  diagram,  the 
locations  of  plastic  zone  boundary  are  plotted  at  different  load  levels,  the  label  on  each 
line  represents  the  level  of  the  applied  traction  in  MPa,  The  hatched  side  of  the  line 
represents  matrix  that  is  currently  plastically  deforming ,  and  we  refer  this  enclosure  as 
the  plastic  zone.  Prior  to  external  loading,  the  matrix  has  partially  plastically  deformed 
around  the  whisker  due  to  plastic  relaxation  of  the  thermal  stresses.  When  a  tensile 
load  is  imposed,  the  shape  of  the  plastic  zone  changes.  It  expands  along  the  longitudinal 
whisker-matrix  interface,  while  unloading  from  the  initial  residual  compression  occurs  at 
the  tip  of  the  whisker. 

As  the  applied  load  increases,  plastic  zones  around  two  adjacent  whiskers  attract  each 
other.  This  is  characterized  by  the  bowing-out  of  the  plastic  zone  boundary  shown  in 
Fig.  2  when  applied  stress  reaches  35.5  MPa. 

When  the  external  load  reaches  about  36  MPa,  the  plastic  zone  boundaries  intercon¬ 
nect  along  BC.  Upon  further  load  increases  ( e.g .  to  44  MPa  in  Fig.  2),  the  only  part  of 
the  matrix  behaving  elastically  is  at  the  tip  of  the  whisker  where  it  is  still  going  through 
unloading  and  loading  in  reverse  direction  with  respect  to  the  compressive  longitudinal 
thermal  residual  stress.  When  the  internal  stresses  gradually  reach  the  matrix  flow  stress, 
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Figure  1.  Schematics  of  an  unit  cell  taken 
from  an  infinite  hexagonally  distributed  ar¬ 
ray.  In  this  unit  cell,  the  boundary  between 
the  initial  compressive  and  tensile  zones  is 
shown. 


Figure  2.  Sequence  of  the  plastic  zone  ex¬ 
pansion  in  the  matrix  of  a  20V%  whisker  re¬ 
inforced  SiC-6061  A1  composite  with  a  ther¬ 
mal  history  of  AT  =  30°  C.  An  external 
load  in  MPa  is  applied  incrementally. 


the  plastic  zone  starts  to  propagate  into  the  initial  compressive  zone.  The  emergence  of 
the  load-induced  plastic  flow  in  the  initial  compressive  zone  then  breaks  the  remaining 
elastic  region  into  two  “elastic  pockets”  —  one  immediately  next  to  the  matrix-whisker 
interface  at  the  tip  of  the  whisker,  the  other  some  distance  away  from  the  corner  of  the 
whisker. 

To  summarize,  the  entire  process  of  plastic  zone  expansion  can  be  divided  into  four 
stages:  (1)  plastic  zone  expansion  in  the  initial  tensile  zone;  (2)  plastic  zone  intercon¬ 
nection  along  the  cell  boundary  BC\  (3)  overcoming  the  compressive  thermal  residual 
stresses  at  the  tip  of  the  whisker;  (4)  fragmentation  of  the  elastic  matrix  to  form  remnant 
elastic  pockets  at  SiCw  tip  surrounded  by  the  matrix  plastic  flow. 

The  matrix  microplasticity  is  also  correlated  to  the  composite  mechanical  properties. 
Fig.  3  shows  the  initial  portion  of  the  stress-strain  curve  produced  by  FEM,  and  how  it  is 
related  to  the  plastic  zone  expansion  process.  At  the  start,  the  stress-strain  relationship  is 
approximately  linear.  However,  because  of  the  thermally  induced  plasticity  in  the  matrix, 
the  thermally  induced  plastic  zone  expands  even  within  this  approximately  linear  region 
as  shown  in  Fig.  2.  As  the  plastic  zones  attract  to  each  other  and  interconnect,  the 
stress-strain  curve  goes  quickly  out  of  the  approximate  linearity.  Once  the  entire  initial 
tensile  zone  is  plastically  deforming,  the  slope  of  the  stress-strain  curve  approaches  a 
constant  again  (Fig.  3).  This  is  due  to  the  fact  that  the  rate  of  propagation  of  the 
plastic  zone  boundary  is  greatly  reduced  because  of  the  ongoing  process  of  overcoming 
the  compressive  residual  stress  in  the  initial  compressive  zone. 

As  the  load  further  increases,  onset  of  the  composite  global  yielding  is  characterized 
by  the  formation  of  the  remnant  elastic  pockets.  At  .2%  yield  stress,  the  entire  matrix  is 
plastically  deforming  except  within  the  elastic  pocket  away  from  the  corner  of  the  whisker. 
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Figure  3.  Initial  response  of  a  20V%  com¬ 
posite  with  AT  =  30°C,  where  the  devia¬ 
tion  from  the  approximate  linearity  is  char¬ 
acterized  by  the  plastic  zone  interconnec¬ 
tion  in  the  initial  tensile  zone. 


Figure  5.  Process  of  the  plastic  zone  ex¬ 
pansion  and  interconnection  in  the  matrix 
of  a  20V%  SiCu,/Al  composites  with  a  ther¬ 
mal  history  of  AT  =  480°C.  The  applied 
stress(MPa)  is  applied  incrementally. 
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Figure  4.  Four  stages  of  the  plastic  zone  ex¬ 
pansion  and  their  chara  cteristics  reflected 
on  the  tensile  behavior  of  a  20V%  SiC^/Al 
composites  with  a  thermal  history  of  AT  = 
30°  C,  a:  expansion  of  the  initial  tensile 
plastic  zone;  b:  interconnection  of  the  ini¬ 
tial  tensile  plastic  zones;  a  deformation- 
induced  plastic  flow  in  the  initial  compres¬ 
sive  plastic  zone;  d:  fragmentation  of  the 
the  remnant  elastic  matrix. 
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Figure  6.  Different  stages  of  plastic  defor¬ 
mation  and  the  tensile  behavior  of  a  20V% 
SiC^/Al  composites  with  a  thermal  history 
of  AT  =  480°C,  a:  steady  plastic  zone  size; 
b:  consumption  of  elastic  matrix  at  the  tip 
of  the  whisker  by  incremental  plastic  flow; 
c:  fragmentation  of  the  remnant  elastic  ma¬ 
trix. 
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Figure  7.  Stress-strain  curves  under  tension 
and  compression  predicted  by  the  FEM, 
where  the  apparent  young’s  modulus  of  the 
composite  is  higher  when  it  is  uder  tension, 
whereas  the  yield  strength  is  greater  when 
it  is  under  compression. 
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Figure  8.  FEM  prediction  of  the 
Bauschinger  stress  factor,  at,  as  a  function 
of  the  applied  forward  strain,  where  o\,  is 
larger  if  the  forward  load  is  in  compression, 
and  the  initial  value  of  <7k  is  negative. 


As  will  be  seen  during  the  remainder  of  this  study,  the  final  location  of  the  remaining 
elastic  zone  is  a  strong  function  of  reinforcement  concentration.  The  correlation  of  each 
stage  of  the  plastic  zone  expansion  with  the  stress-strain  relationship  is  shown  in  Fig.  4. 

For  AT  =  480°C,  which  is  a  more  realistic  temperature  change  when  cooling  from  the 
annealing  temperature,  the  evolution  of  the  matrix  plasticity  is  shown  in  Fig.  5,  where 
the  same  convention,  as  employed  in  Fig.  2,  is  used.  Due  to  a  larger  AT,  the  thermally 
induced  plastic  flow  has  spread  over  the  entire  matrix  prior  to  any  external  loading.  As 
compared  with  the  cooling  of  AT  =  30°C,  the  development  of  matrix  plasticity  due  to  a 
tensile  load  is  the  same  in  characteristics  except  that  the  first  two  stages  of  the  plastic 
zone  expansion  occur  during  “thermal  loading”,  which  leads  to  that  the  changes  of  the 
matrix  plasticity  is  concentrated  at  the  tip  of  the  whiskers. 

Fig.  6  shows  the  predicted  stress-strain  curve  of  the  composite  along  with  the  corre¬ 
sponding  characteristics  of  the  plastic  zone.  Comparing  with  Fig.  3,  the  “proportional 
limit”  is  considerably  smaller,  the  difference  comes  from  the  fact  that,  for  AT  =  480° C, 
the  matrix  is  experiencing  a  different  stage  of  plastic  zone  expansion.  When  the  initial 
loading  is  applied  within  the  approximate  proportional  limit,  the  matrix  is  unloading 
in  the  initial  compressive  zone,  where  the  boundary  of  the  plastic  zone  is  insensitive  to 
the  external  deformation.  However,  it  takes  significantly  more  additional  work  to  spread 
the  deformation-induced  plasticity  to  reach  the  stage  of  fragmentation  of  remnant  elastic 
matrix.  This  is  due  to  the  increase  of  compressive  residual  stress  and  the  amount  of 
thermally  induced  work  hardening  when  AT  is  larger,  i.e.,  the  thermally  induced  dislo¬ 
cations  due  to  ACTE. 

3.2  COMPOSITE  ASYMMETRIC  CONSTITUTIVE  BEHAVIOR  AND  THE  THER¬ 
MALLY  INDUCED  PLASTIC  FLOW 

Fig.  7  displays  the  uniaxial  stress-strain  curves  generated  by  FEM  when  a  20V%  com¬ 
posite  is  under  tensile  and  compressive  loads,  respectively.  Fig.  8  shows  the  Bauschinger 
stress  factor,  <76,  as  a  function  of  the  applied  forward  strain.  Here,  the  Bauschinger  stress 
factor  is  defined  as  <?6  =  aF  —  oR,  where  aF  and  aR  represent  the  forward  and  reverse 
flow  stresses,  respectively.  In  this  figure,  <76  is  consistently  larger  when  forward  applied 
stress  is  in  compression.  These  results  match  the  observations  by  Arsenault  et  al.  [3]. 


Figure  9.  Matrix  plastic  zone  expansion  Figure  10.  Schematics  of  the  morpholo- 
process  of  a  20V%  SiC^/Al  composite  un-  gies  of  the  plastic  zone  under  different  load- 
der  compressive  loading  condition  with  a  ing  conditions,  where  tensile  loading  will 
thermal  history  of  AT  =  480° C.  produce  a  pseudo-isostrain  condition  and 

compressive  loading  will  produce  a  pseudo¬ 
isostress  condition. 


In  addition,  FEM  analysis  also  predicts  that  the  Bauschinger  stress  factor  is  initially 
negative  (a b  <  0),  which  means  that  there  is  an  initial  reverse  strain  hardening  rather 
than  normally  observed  reverse  strain  softening. 

3.2.1  Asymmetric  Apparent  Young’s  Modulus 

It  is  commonly  believed  that  the  asymmetric  composite  apparent  Young’s  modulus  is 
related  to  the  asymmetric  development  of  the  shear  stresses  at  the  interface  near  the  tip 
of  the  whisker  [11,  12].  From  FEM  modeling,  Levy  et  al.  found  that  the  interfacial  ther¬ 
mal  residual  shear  stress  near  the  end  of  the  whisker  increases  with  compressive  loading, 
and  decreases  with  tensile  loading  [12].  Therefore,  they  argued  that  micro-plastic  flow 
induced  by  the  asymmetric  development  of  the  shear  stress  near  the  longitudinal  interface 
produces  early  micro-yielding  during  compressive  loading,  so  that  the  apparent  Young’s 
modulus  is  lower.  However,  this  theory  cannot  explain  how  the  asymmetric  development 
of  matrix  shear  stress  leads  to  a  higher  compressive  composite  yield  strength  as  shown 
in  Fig.  7. 

On  initial  tensile  deformation,  the  development  of  the  interfacial  shear  stresses  induces 
unloading  in  the  matrix.  However,  from  the  evolution  of  the  matrix  plastic  flow  shown 
in  Fig.  5,  the  effect  of  such  unloading  is  restricted  to  the  immediate  proximity  along 
the  longitudinal  interface  near  the  tip  of  the  whisker,  and  its  influence  remains  at  small 
applied  load,  e.g.,  at  a  load  as  small  as  2.2  MPa,  the  plastic  zone  starts  t*<  expand  toward 
the  tip  of  the  whisker. 

Fig.  9  shows  the  characteristics  of  the  plastic  zone  under  a  compressive  load.  Compar¬ 
ing  the  morphologies  of  the  plastic  zone  under  tensile  and  compressive  loading,  respec¬ 
tively,  as  shown  in  Figs.  5  and  10,  the  unit  cell  may  be  divided  into  three  distinct  regions: 


elastic,  plastic  deforming  regions  and  the  whisker.  Fig.  10  display  an  idealization  based 
on  this  partition.  If  we  consider  the  two  stiffer  regions  to  be  an  imaginary  phase,  we 
may  define  the  morphology  of  the  initial  development  of  matrix  plastic  zone  under  ten¬ 
sion  and  compression  as  pseudo-isostrain  and  pseudo-isostress  conditions,  respectively, 
in  reference  to  the  isostrain  and  isostress  loading. 

For  isostrain  condition: 

Ec  =  E,(l  -  Vj)  +  EjVj.  (1) 


For  isostress  condition: 


Ec  = 


_ Ei  Ei _ 

(1  -V3)E2  +  ElV3‘ 


where  £,  (i  =  1,2)  and  Ee  are  the  Young’s  moduli  of  constituent  phase  t  and  the 
composite,  respectively.  Vi  is  the  volume  fraction  of  constituent  phase  t. 

For  the  pseudo-isostrain  model  (Fig.  10),  the  stiffness  can  be  evaluated  by  substituting 
Eq.  1  with  the  equations  shown  below: 

p  _  _ EwEm _ 

1  V;E»  +  (1-  V°)Em'  w 


where  V°  corresponds  to  the  volume  fraction  of  the  plastic  zone  in  the  pseudo-isostress 
model;  Em  and  Ew  are  the  Young’s  moduli  of  the  matrix  and  whisker,  respectively;  Emp  is 
the  matrix  work  hardening  rate.  Substitute  Eqs.  3  and  4  into  Eq.  1  and  consider  V2  ■=  V* 
where  V*  corresponds  to  the  volume  fraction  of  the  plastic  zone  in  the  pseudo-isostrain 
model,  the  modulus,  E ',  of  the  pseudo- isostrain  model  can  be  obtained: 


El  = 


(1  -  V‘)EwEm 
V;E„  +  (l-V°)Ev 


?  uf 

Jmp  v  p 


For  the  pseudo-isostress  model  (Fig.  10),  Ei  =  (1  -  V*)EW  +  V‘Em,  E2  =  Emv,  and 
V2  =  V° .  Substitute  E\  and  E2  into  Eq.  2,  we  obtain: 

„  =  U-V:)E.Emr  +  V;EmEmr 
'  (1  -  Vf)Emp  +  (1  -  Vf)V;B.  +  V;V;E„  u 

where  E°  is  the  modulus  for  the  pseudo-isostress  condition.  Considering  material  prop¬ 
erties:  Em  =  68.3  GPa;  Emp  =  2.08  GPa;  Ew  =  483  GPa  [9],  and  for  a  unit  cell  shown  in 
Fig.  1  which  corresponds  to  a  20V%  composite,  we  obtain  £f  =  9.887  GPa;  El  =  50.09 
GPa,  i.e.,  Elc  >  E°.  The  magnitude  of  the  asymmetry  of  the  Young’s  modulus  is  ex¬ 
aggerated  because,  as  compared  with  that  of  the  actual  plastic  zone  shown  in  Fig.  9, 
the  morphology  of  the  plastic  zone  shown  in  Fig.  10  during  compression  leads  to  a  lower 
bound  for  the  composite  stiffness. 

From  this  result,  it  can  be  seen  that  the  asymmetric  apparent  composite  Young’s 
modulus  is  due  to  an  asymmetric  plastic  zone  expansion  which  induces  isostress  or  isos¬ 
train  type  of  loading.  This  process  continues  until  the  morphology  of  the  plastic  zone  can 
no  longer  be  represented  by  the  idealization  shown  in  Fig.  10,  that  is,  when  the  applied 
stress  is  approaching  the  composite  yield  stress. 


3.2.2  Asymmetric  Bauschinger  Effect 

In  explaining  the  phenomenon  of  asymmetric  Bauschinger  stress  factor,  the  “back 


stress”  model  can  be  modified  [13],  in  which  the  forward  (of)  and  reverse  (erf)  flow 
stresses  can  be  expressed  as: 

\af  |  =  oy  4-  A  erf  +  |A<X/fc|  (") 

|<xf  |  =  (Ty  -  Atrf  +  |A<7/*|  (8) 

where  oy  is  the  yield  strength;  A  erf  and  A  erf  are  the  contributions  of  the  back  stress 

during  forward  and  reverse  cycle,  respectively;  A <T]h  is  the  plastic  friction. 

From  Eq.  7  and  8,  the  Bauschinger  stress  factor  can  be  expressed  as:  o\,  =  \of\  — 
|  erf  |  =  A  erf  +  Aerf .  Approximating  the  back  stress  as  the  matrix  residual  stress,  we 
obtain  the  following  equations: 


where  ( aTm )  is  the  thermal  residual  stress;  C(eF)  is  a  function  of  forward  applied  strain, 
eF .  ( a represents  the  changes  of  residual  stress  induced  by  the  external  deformation. 
Then  the  Bauschinger  stress  factor  is 

<7»«-[l+C(EF)]  j^jto  (11) 

Considering  C(eF)  =  1  -  keF  +  0(eF),  where  k  >  0,  and  0(eF )  represents  the  terms 
with  higher  order  than  eF.  For  the  first  order  approximation,  0(eF)  may  be  neglected, 
that  is,  C(eF)  =  1  —  keF .  Substituting  C  into  Eq.  11,  we  obtain: 

^  =  (^-2)^(0  (12) 


Eq.  12  reflects  the  general  trend  of  the  changes  predicted  by  the  FEM  shown  in  Fig.  8. 
For  tension-compression  loading  scheme  (er  >  0), 


Ob 


( 


<  0 
>  0 


eF  <  1- 

£F  >  - 
t  —  k- 


(13) 


For  compression- tension  loading  scheme  (eF  <  0),  ex*,  >  0.  In  both  cases, 


k  >  0. 


4  CONCLUSIONS 

From  the  above  analysis,  the  following  conclusions  can  be  reached: 

(1)  The  evolution  of  the  matrix  plasticity  in  a  whisker  reinforced  SiC/Al  composite 
is  characterized  by  the  expansion  and  interconnection  of  the  thermally  induced  plastic 
zones. 

(2)  The  process  of  plastic  zone  expansion  can  be  divided  into  four  different  stages. 
Plastic  zone  expansion  in  the  initial  tensile  zone;  plastic  zone  interconnection  in  the 
initial  compressive  zone;  overcoming  the  longitudinal  compressive  matrix  thermal  residual 
stress  at  the  tip  of  the  whisker;  and  the  fragmentation  of  the  remnant  elastic  matrix  into 
pockets. 


(3)  The  onset  of  global  yielding  is  characterized  by  the  fragmentation  of  the  remnant 
elastic  matrix  at  the  tip  of  the  whisker,  where  isolated  remnant  elastic  pockets  remain. 

(4)  The  composite  asymmetric  constitutive  behavior,  the  asymmetric  Bauschinger  ef¬ 
fect  and  the  apparent  Young’s  modulus,  can  be  related  to  the  asymmetric  matrix  plastic 
flow  due  to  the  thermally  induced  plasticity  and  residual  stresses. 
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